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ABSTRACT 



The reliaDllity of the structure can be calculated by combining the 
non-uniform stress distribution of a structure from finite element analysis 
with the statistical strength theory under two dimensional non- 
homogeneous, uniaxial stress limitations. The specimen models with and 
without notch are the sample structures to illustrate the calculations. The 
statistical strength of the structure is cast in the standard Weibull form 

characterized by the structural scale parameter, and shape parameter, 
cx^ which are functions of the material scale and shape parameters (<x and 

3.) Therefore, this thesis demonstrated that the scale parameter of the 
structure remains constant for different load magnitude, only the 
structural geometry and loading condition cause the change of material 

scale parameter, ^ to the structural scale parameter, 3^- The results also 

shows that the degree of uniformity of stress within structural element 
and the linearity of the load magnitude and stress distribution within the 
structure affect the accuracy for the calculation 
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I. INTRODUCTION 



A. RELIABILITY REQUIREMENT FOR COMPOSITE STRUCTURES 

R 2 liability theory became an independent scientific discipline at the 
beginning of the 1950's under the influence of the rapid development of 
radio electronics, computer, devices, and rocketry. Modern radio electronic 
apparatus and digital mechanics consist of a very large quantity of 
components. If the failure of one component leads to the failure of the 
unit as whole, then evidently the possibility of high reliability (freedom 
from failures) will diminish rapidly as the unit becomes more complex. In 
this connection, problems arise on the prediction of the reliability of 
planned units, on the development of measures to increase the reliability, 
on a fundamental for reliability testing methods, etc. All these questions 
are the subject of reliability theory. 

The mechanical strength and stiffness of a structure are one of the 
aspects of reliability. The engineering design concept is not only solving 
the stresses and strains which originate in structures under various 
external imposed conditions but also trying to make them to be operated 
sufficiently reliably throughout their established utilization time. 
Therefore, in the concluding stages of an engineering design structural 
mechanics Inevitably comes into resolution with reliability theory. For 
critical application of composite as an aircraft structure, the study of 
reliability in composite structures is required. 
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B. HISTORICAL BACKGROUND 



Bullock [ Ref. 1 1 used Weibull theory to predict the strength of the 
composite materials between different laboratory sample configurations. 
He investigated the strength ratios of composite materials in flexure and in 
tension. He assumed the shape parameter (the exponent characterize 
flaw-density that determines the scatter of strength of the material) was 
not changed when the geometry or size of the material was changed. He 
used the same shape parameter for the tensile test and flexure test which 
based on homogeneous state of stress and heterogeneous state of stress 
respectively. He found, however, based on the experimental results, the 
hypothesis that the shape parameter was constant was not realistic. 
Changes in shape parameter were in fact observed. Whitney and Knight, n 
I Ref. 2 1 used a statistical strength theory based on Weibull distribution to 
explain the difference between unidirectional tensile data generated from a 
flexure test and a standard tensile coupon. The result was shown a 
significantly larger variation in tensile strength versus flexural strength! 

Rosen ( Ref. 3 1 presented a theoretical and experimental treatment of 
the failure of composites, consisting of a matrix stiffened by uniaxial ly 
oriented fibers, when subjected to a uniaxial tensile load parallel to the 
fiber direction. He observed that a portion of the fiber at each end was not 
fully effective in resisting the applied load because of the axial load was 
transmitted by shear through the matrix to adjacent fibers. Basically, the 
fibers failed as a result of a statistically distributed flaws or 
imperfections and composites failed as a result of a statistical 
accumulation of such flaws over a given region. Therefore, the ineffective 
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length of the fiber was introduced in order to consider the composite to be 
composed of a series of layers of this length. Then the segment of a fiber 
within a layer might be considered as a link in the chain that contribute the 
fiber. Each layer was then a bundle of such links, and the composite as 
whole can be modeled by series of such bundles. Weibull theory was used 
to define a statistical distribution of flaw or imperfections that result in 
fiber failure under applied stress. The statistical accumulation of such 
flaws within a composite material was demonstrated to be cause of 
composites failure. He concluded that the stress concentrations occurred 
in fiber adjacent to the fibers break were considered only to balance the 
effect of the variation of ineffective length which affect to the composite 
strength. Knight, C.E. [ Ref. 4] used a finite element stress analysis and 
statistical strength theory to assess the Influence of the stress 
concentration on the ultimate strength calculation for composites from 
the Split-D Test. He observed that the geometry changed causes the 
nonlinearity and the boundary condition would be changed. 

The Probabilistic Statistic Failure of Composite Materials [ Ref. 5 ] 
demonstrated how to develop a theory that combined the statistics of 
composite material failure with the orthotropic nature of composite 
materials. The effects of loading history and the probabilistic location of 
the failure for a composite can be accounted using this theory. In order to 
improve the efficiency of composite designs for structural analysis, a 
theory that combining the statistics of the experimental failure data with 
an orthotropic material failure response was developed. The analytical 
theory developed extent deterministic failure laws to include the 
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statistical scatter observed experimentally. Tsai and Wu [ Ref. 6 ] are one 
of the investigators who described a general theory for the failure of 
composite materials as the strength tensor theory but determining is 
complicated by wide range in failure properties of composite materials and 
scatter in experimental results. The Weibull weakest link , theory was 
developed and described a description of the spatial distribution of the 
failure locations with a sample calculation performed for the case of a 
uniformly loaded tensile specimen. Three point bending specimens were 
used to predict failure location by applying the basic weakest link theory. 
Four point bending specimens were used to illustrate the effect of loading 
history which were evaluated by applying numerical integration to the 
probability of failure density distribution curves. The compressive failure 
load was demonstrated that it is not deterministic failure and also having 
some scatter which is independent of tensile scatter so that the Weibull 
failure distribution for tension and compression are different. 

In ordinary structures having many stress components active, a 
decision on the probability for survival must be based on statistical 
descriptions that admit multiaxial stress. The development of combining 
weakest link theory with the strength tensor representation is required. 

C. OBJECTIVE AND SCOPE 

The driving motivation for this research is that many critical 
components of modern aircraft structure are made of composite materials. 
However, the reliability design is necessary for engineering designs 
requires not only the stress analysis but also statistical theory of failure 
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strength of composites, The latter motivated the investigation of 
statistical theory for composites failure strength where are summarized in 
the last sectioa Any aircraft structural components inevitability have the 
non-uniform stress in the neighborhood of edge-notches and holes. Size 
effect is the one that is importance for the composite study because of the 
difficulty and the economy that one could not experiment the several large 
samples to predict the reliability of the components. Even the 
experimental test samples can be made, the difficulty to identify the 
failure location is also impossible. Finite element method and the 
statistical theory can be combined to explore the possibility of prediction 
both the reliability of the structure and the failure location. Basically, 
finite element method can be used to calculate the stress distribution that 
occurs in the structure under an arbitrary loading boundary condition. The 
structure is divided into elements be sufficiently small that the stress is 
uniform. This solution of physically continuous system can then be 
replaced by the solution of the discrete system of the elements. The 
statistic Weibull theory can be used to model composite materials 
composing of a series of bundles which are in turn composited series of 
fiber links. Another assumption for the statistical Weibull theory is the 
local stress in each link must be uniform. For the connection to the two 
dimensional finite element method, instead of using link, the element will 
be introduced. For one dimensional case, the link can be used but in two 
dimensional case, the element has a width. The operation of size effect 
formulation based on Weibull theory and finite element method will be 
presented in next chapter. 
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The scope of this research limits to the treatment of statistical 
strength of composite structures which are two dimensional, subjected to 
nonhomogeneous uniaxial stress including size effect. The specimen 
models are plate with and without notch subjected to the tensile 
displacement loading case. The finite element code name NIKE2D [ Ref. 7 ] 
is used. NIKE2D is a two dimensional finite element program which is 
developed by Hallquist, J.O. [ Ref. 7 ]. It can be used to calculate linear and 
non-linear finite element problems. In this thesis, only the linear portion 
is utilized. MAZE I REf. 8 ] computer program is used to generate specimen 
model, grid, define boundary conditions, and also create the input file that 
can be immediately use for NIKE2D program. ORION [Ref. 91 is the post 
processing of NIKE2D which can be used to create contour curves of stress, 
displacement, etc. The observation of uniform stress in each element will 
be made. The observation of using half of the entire model by symmetrical 
will be made, which bases on the idea that the model can be divided into 
smaller elements, when compare to the entire model. The linear 

relationship between external loads and the local stress in each element 
will be verified. The plate without notch is used as bench mark to compare 
the reliability results to the plate with notch. The specimen model 
dimension bases on specimen dimension which using in the shear test 
experiment. The shear model that involves with the multiaxial stress 
which can not yet be treated in this thesis. Future work requires 
generalization of current results to multiaxial stress. When such extension 
is made, realistic prediction of composite structures will then be possible. 
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II. STATISTICAL MODEL OPERATION 



A. REVIEW OF WEIBULL WEAKEST LINK FAILURE THEORY 

consider an infinitesimal volume element SV at point x and subjected 

to stress state Ojj. The probability of failure Sf of volume element 8V 

is assumed to be given by: 

Sf = >? [ [x] ] SV ( 2.1 ) 

where: 'P = failure parameter which is a function of Ojj which is a 

function of position, x. 

It is assumed there are only two possible states for the volume 
element, SV , either the volume element, survives the loads or fails. From 
this assumption the probability for survival or reliability is simply: 

SR = 1 - Sf = 1 - 'P[ OjjMlSV (2.2) 

The basic assumption in weakest link theory is: if any volume, SV , in 
the total structure volume, V , fails, then the entire structure occupying 
volume V fails. 

Consider the structural volume, V, to be divided into N volume 
elements SV|^ , with each of the volume elements located at point X|^ . 

Then the total volume is given by: 

V= lim S SV,^ = Jy dV: k = I, 2 . . . N (2.3) 

where: n ♦ 
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The probability for the volume V to survive, R , under weakest link 
theory, is given by: 

R = lim n SR,^ = lim n { 1 - ] SV,^ } ( 2.4 ) 

where: k = 1, 2 . . . N , and 

N-* «> . 

Equation 2.4 can be equivalently written as: 

lnR= lim Z In { I - 'J'f Ojj [X|^] SV|^ } (2.5) 

where: k = 1, 2 . . . N , and 

N- 00 

The quantity is small and from the Taylor series expansion 

In (1 + x) ^ X for X « 1 

Equation 2.5 becomes: 

lnR = lim Z -{ 'J'l Ojj[X|^]]SV|^} (2.6) 

or R = exp -{ Jv n <Tjj [x,^] ] dV } ( 2.7 ) 

where; k = I, 2 . . . N , and 
N -*00 

Equation 2.7 is the basic equation derived under the assumption of 

weakest link theory. Note that the stress distribution need not to be 
uniform; the material can be generally anisotropic and the material can be 
non-homogeneous. The equation is limited in that it describes only the 
reliability of the entire structure and provides only a quantitative 
description of where the structure will fail. The reliability becomes: 

R = exp - { Ojj/3 }‘^ 
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where: ex = shape parameter, and 

3 = scale parameter 



In the next section, basic theory will be used to obtain the 
formulation of size-effect for one dimensional, non-homogeneous, 
uniaxial stress. 

B. SIZE-EFFECT FOR NON-HOMOGENEOUS STRESS. ONE DIMENSIONAL, 
UNIAXIAL STRESS 

The value of (X|_j and 3 li are estimated by the experimental 

procedure from the given n numbers of sample length LI . For a sample 

length L , the evaluation of and 3[_ can be implemented according to 
the following procedure. 

There are j links in length L (Figure 2.1-a) each links are not 
necessarily of equal length, but the stress in each link is homogeneous, ie 
the stress cfj^P within the p^^ link satisfies the following condition: 



xa^P-(/|p xa^dxl/lp ? t 



( 2 . 8 ) 



Where: 



€ can be any arbitrarily small value 
Ip = length of p^^ link 



P = 1. 2 . . . j 
The reliability for each link is: 



Rp(o/) = I - f(a/) = exp -(o//3p)«p 



(2.9) 
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Figure 2.1 



(a) One Dimensional Links Configuration 

(b) Two Dimensional Specimen Configuration 
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There are (lp/l) unit links within Ip , the reliability of each unit link 
within Ip is: 

Rp(OjjP) = exp ( 2.10 ) 

where: « and ^ for unit link are material constants 

The reliability for each link Ip is: 

Rp(Ox^) = IT {exp ; m = I. 2 . . . (Ip/1) 

= exp {-(lp/l)(o 

where: (ip/l) is the number of unit link in Ip 

The reliability for then entire length L is: 

Rl " R] (^ 1 ) ^ 3 (^ 3 ) R 4 (<T^) . . . Rj(<Tj) 

= exp {-(li/1)(<3 exp (-(l 2 / 1 )(o 2 / 3 )^. . . exp {-(lj/1)(o 
= exp -{(li/1)(o 1 / 3 )^ + (l 2 /l)(o 2 / 3 )*^ + . . . + (lj/1)(o j/3)*^ 

Rl = exp -{ S [(lp/1)(0 //3)‘^ } : p = 1, 2 . . . j ( 2.11 ) 

Under uniform stress condition within an element, let Op is the 

element stress which is constant within Ip^^ length element and Op is 

calculable from a given boundary condition (B.C.) characterizable by (P), 
(eg, a load vector). The stress due to the same B.C. with a different 
magnitude P Is: 

Op[P] = Op[{P}] (P/{P}) ( 2.12 ) 

where: P can be of any arbitrary magnitude in linear system, or 

P is of neighboring magnitude of local linear system. 
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<Jp[{P}l is expressible from structural mechanics calculations, 
(eg, a finite element method.) 

In order to reduce the value of which is a function of P in 
equation 2.11 to a single random variable {P} , the substitution of cfj^P[Pl 

by is made within the same manner as the expression in equation 

2.12. Thus, the equation 2.11 becomes: 

Rl = exp -{ Z l(lp/l)(a /[{P}l/{P})'^P/3)‘^] } ; p = I. 2 . . . j 

Expressing the above equation into the Standard Weibull form: 

R|_ = exp -{ (P/3i_)®l } 

Then: = « 

= ?/ ( s: [(lp/l)(o /1(P)1/{P))“) ; p = I, 2 . . . j 

Let: F = ( 2 l(lp/l) (cr/[(P}l/{P))l“ ; p = I, 2 . . . j 

Thus: Sl = ?F 

That is, given a structure of length L subjected to boundary condition 
parameter {P} ( which may be of a single or multiple load Pi , P 2 , . . . 
vectors characterizable by a single parameter {P} ) which give rise to a 
unidirectional stress Op computable from the methods of mechanics ( eg, 

finite element analysis ). Ip are element lengths which are segmented as 
small as necessary to assume the stress to be uniform. Upon sampling n 
structures measuring failure loads. 
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oc|_ , 3j_ are shape and scale parameters ( in dimension of load 

estimated from D(Pp^) and structural size L . 

o< , 3 are shape and scale parameters ( in dimension of stress which 
are geometric independent material parameters ). 

F is a conversion factor calculable for structural ( finite element ) 
analysis. 

For special case, homogeneous, one dimensional uniaxial stress: 



F = ( E l(ip/i) (a/l(P)l/(P)))°' 


; P = 1. 2 . . . j 


'* = 0 and (P) = Pn: 




F = ( (a/Pn)« { E (Ip/I) '''« }« 


; P = 1. 2 . . . j 


= 1 (O/Pp) { E (Ip/I) 


; P = 1, 2 . . . j 


= (0/Pp)-'(L/l)-'''« 





Thus: 

R|_ = exp -{ [ (P/3)/(a/Pn)-' /(L/l) )<x ) 

= exp -{ (L/l) 1 (C/3)(P/Pn) 1“ ) 

This recovers the well known classical result demonstrates the 
consistency of the current formulation is true. 

For each given non-homogeneous boundary condition, a corresponding 
correction factor F can be calculated by finite element such that the size 
effect can be predicted. By weakest link in length X ; it is inferred that 
the plane of facture is perpendicular to X . That means the failure plane 
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has no thickness. Since the one dimension specimen has no width, fracture 
width is not a parameter. 

C. SIZE-EFFECT FOR NON-HOnOGENEOUS STRESS, TWO DIMENSIONAL, 
UNIAXIAL STRESS 

This is the extension of one dimensional, uniaxial stress to two 
dimensional uniaxial stress. The definition that uniaxial stress implies 
failure has zero dimension in X direction is adopted for one dimensional, 
uniaxial stress , which has been described in the previous section. For two 
dimensional, uniaxial stress, failure automatically requires definition of 
the failure plane because the specimen has dimension in Y direction 
(Figure 2.1-b), and it is no longer conceptually a point as in the one 
dimensional case. In the similarly procedure as in one dimensional case, 
the specimen is divided into sufficiently small areas within which the 
stress in direction of force is assumed to be uniform. It is seen that not 
only the length of the specimen is divided into small element lengths but 
also the width must be divided into small element widths in order to 
provide element areas which are satisfies uniform stress assumption. Thus: 

for all Ip 

Where: p = I, 2 . . . j 

a°x'^ ■ wq/ yOx^a^Wq « fy 3" «q 

Where: q =1, 2 . . . k 

Weakest link in X direction is defined by Eq. 2.8 and unit length in X 
direction is defined by Eq. 2.9. 
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Weakest link in Y direction is: 

where: Oj^P^ is stress in the q^*^ element, and q = I, 2 . . . k 

This means that If any one of the q^^ fall, all elements fail. However, 
the unit (metric) width can not be arbitrarily defined. It is visualized that 
within each metric area, there Is one domineering crack (largest crack 
perpendicular to Oj^). When subjected to the domineering crack 

extends without bound. From fracture mechanics: 
kc=a/a7 

where: k^ Is a material constant and can be independently measured. 

a^. Is weakest link in Y direction, which can be used as the unit 
metric width. 

Thus: a^P^ = (k^Oj^P^)^ for element pq*^ 

As well as the division of length element Ip into unit length, the 

division of width element is necessary. In stead of unit length 1 , a^P^ is 

used in order to indicate the critical width in each element, thus, the 
reliability of entire element is: 

f’pq'®/'’' = ■< Cpq/'Xw pq/3(.P''Xa/<l/3)« ) 

= Bxp -( (lpq/l)lw pq/(k^a^P'l)2l(a)jP<t/{)« } 

= exp -{ (Ipq/lXw pq/kp2)(a,P<l)2*«/J« } 
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where: pq = 1, 2 . . . jk ( subscript is changed from two dimensional 

array into one dimensional array.) 

The reliability of entire specimen (which is the structure) is: 

■ R| R 3 ■ • • Rjk 

= exp -( Z (Ipq/lXw pq/K2)(aKP1)2*«</J« ) , pq = I, 2 . . . jk 
From the linear relationship between external load and local stresses: 
o/^[P] = o/^[{P}] (P/{P)) 

where: a^P^[{P}] is element stress which is calculated under certain 

loading boundary condition while the magnitude of load is {P}. 

a/^[Pl is element stress which subjected to the same 
boundary condition as above but different in magnitude, the magnitude is P. 

Thus: = exp -{ E (lpq/l)(w pq/kc2)((P/{P})a/^[{P}])2^«/^« } 

where: pq = 1, 2 . . . jk 

Express it in term of the Standard Weibull for the structure, thus; 

Re = exp -{ (P/^e)"^ ( 2.13 ) 

where: o<^ = <x + 2 

{ Z (Ipq/lXw pq/k(.2Xa/'’((P)l)‘>'*2 

where: pq = I, 2 . . . jk ( 2.14 ) 

The values of « and 3 of material are altered by the boundary condition 
and the geometry of the specimen (or structure) changed but they do not 
depend on the changed of the magnitude of the load. Therefore, the shape 
parameter, <x and the scale parameter, 3 can be solved for any given 
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structural geometry and boundary conditioa The finite element method can 
be used to calculate the stress distribution of the specimen by dividing the 
specimen into small elements within which the stress element is uniform. 

For the simple limiting 2-D case, the specimen without notch, the 
stress distribution for entire specimen is uniform. Thus: 

= Ojj = constant 

{ j. )(-l/c<*2) 

= (P) Kc(2/o<e) 

Where: A = E [ IpqWpq ] , pq = I, 2 . . . jk 

(Xg = (x+2 

Thus: Rg = exp -{ (A/k^,^) (p/{p}) ( 2.15 ) 
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III. APPROACH TO THE PROBLEII 



A. SPECiriEN MODEL CONFIGURATION AND COORDINATE SYSTEM 

The specimen model has 1.99 " in length and 0.968" in width is align 
in the y-z plane, (cartesian right hand coordinate system is used). The 
configuration and the coordinate of the specimen without notched and the 
plate with 90 degree notch configuration are shown in Figure 3.1. The 
fiber direction which makes the orthotropic material having large elastic 

modulus ( E^^ ) is considered to be along the y-axis (Figure 3.2). 



B. FINITE ELEMENT PROCEDURE 
1. Material Constants 

The finite element program requires defining the material 
constants [ Ref. 7 ] for orthotropic elastic to be used in stiffness matrix 
calculatioa The material constants of typical graphite epoxy composite 
which used in this thesis are: 

E|j = 18.65 Msi (Mega Pound per square inch) 
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Figure 3.1 
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A = i.99' B= 0-968" 

Dimension of Specimen Models without Notch and with 
Notch 
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Figure 3.2 Specimen Models Coordinate System 
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G ,2 = 0.84 nsi 



2. Boundary Conditions 

The specimen model is subjected to the displacement boundary 
condition (B.C.) in order to be able to use the entire specimen model in the 
calculation. An original problem the specimen model is subjected to the 
displacement B.C. value u^ (Figure 3.3) on the right hand side in the 
positive y-direction and the left hand side is subjected to the displacement 
B.C. value in opposite direction. In order to fix the specimen model 

from rigid body displacement, the left side of the model must be constraint 
in the direction of load. Therefore, the displacement value on the right 
side of the model becomes two time of u^. The transform of the 

displacement B.C. is shown in Figure 3.4. 

3. Data Planning and Purposes 

NIKE2D is two dimensional axis symmetry finite element code 
that the thickness is assumed to be infinite. The plane strain analysis is 
used by assuming no displacement along x-axis and no rotation around y and 
z axis, which It is satisfy for this purpose. There are eight finite element 
program data output to be discussed in this thesis: 

3.1 Specimen model without notch subjected to the 
displacement B.C. value 0.01” 

3.2 Same specimen as 3.1 but the displacement B.C. value is 
O.OOr. The purpose of these two cases is to demonstrate that there is no 
stress concentration along displacement direction and the linear 
relationship between displacement and local stress. 
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Figure 3.3 Displacement Boundary Condition 
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Figure 3.4 Transformation of Displacement Boundary Condition 
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3.3 Specimen model with notch subjected to the displacement 

B. C. value 0.01" and the stress calculation bases on four-point integration 
(Ref. 7 ) in order to observe the uniform stress in each element. 

3.4 The same as 3.3 but the calculation bases on one-point 
integration at the centroid of each element. 

3.5 Half of the specimen model with notch subjected to the 
displacement B.C. value 0.01’ and four-point integration is used. This case 
will be used to compare with the entire specimen model in order to insure 
that the half specimen can be used because of the symmetrical 
configuration. 

3.6 Same as 3.5, but one-point integration is used. 

3.7 Same as 3.6, but the displacement B.C. value is 0.004. 

3.8 Same as 3.6, but the displacement B.C. value is 0.001. 
One-point integration results will be used to be input data for 

statistical program. 

C. STATISTICAL PROCEDURE 

After obtaining the results from the finite element program, 
four-point integration results have to be considered to insure that the 
tensile stress in each element is uniform or almost uniform. If they are 
not uniform, the size of the element need to be reduced. However, due to 
computer memory limitation, absolute uniform stress may not be 
attainable. One should realize that what spatial domain that the element 
size should be reduced and what area the element need not be changed. The 
symmetrical of the model can be considered to be used only half of the 
entire model in order to increase the number of elements. 
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Th2 statistical siz2 2ffect in two <Jim2nsional, non-homog2n20us, 
uniaxial str2ss program in fortran 77 is writt2n to b2 us2d to calculat2 th2 
r2liability of th2 sp2cim2n mod2ls by using th2 finit2 2l2m2nt r2sults and 
shape parameter, o< and scale parameter, 3 from several numbers of certain 
sample size. There is one program including: 1 ) the specimen model 

without notch, 2 ) the entire specimen model with notch, and 3 ) the half 
(top side) specimen model with notch. 

This program is presented In Appendix A. There are 12 output data 
files which are shown in Table 4 . 3 . 
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IV. DISCUSSION AND CONCLUSION 



A. EXAHINATION FOR UNIFORHITY OF UNIFORM STRESS WITHIN 

ELEMENT 

The four-point integration are obtained in entire and half specimen 
models subjected to the displacement boundary condition (B.C.) value 0.01. 
For the full specimen model, the element number 544 (shade in Figure 4.1) 
is considered, the tensile stress values can be read In figure 4.2 ,where the 
stresses of this element are tabulated under rows 2173, 2174, 2175, and 
2176 (relatable to the element by dividing the last number by 4 recovering 
544) and the tensile stresses are 0.03663 Msi, 0.03738 Msi, 0.01933 Msi, 
and 0.0 1851 Msi (Figure 4.2) respectively, which demonstrated that the 
stress within element is not exactly uniform but the differences are small. 
The element number of the half specimen model at the same position as 
element number 544 of the full specimen is 442 (shade in Figure 4.3) and 
the tensile stress can be read in Figure 4.4, where the row numbers are 
1765, 1766, 1767, and 1768 and the tensile stresses are 0.03667 Msi, 
0.03742 Msi, 0.01937 Msi, and 0.01861 Msi respectively, again they are not 
exactly uniform but the differences are small. The element number 192 in 
Figure 4.5 (full specimen model) is outside the notch region, the tensile 
stress values can be read in Figure 4.6, where the row numbers are 765, 
766, 767, and 768 and the tensile stresses are 0.06891 Msi, 0.06888 Msi, 
0.06852 Msi, and 0.06855 Msi respectively , which shown that the stress 
within element is not uniform but the differences are more uniform than 
element in the notch region. The element number of the half specimen 
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model at the same position as element number 192 of the full specimen is 
230 (Figure 4.7) and the tensile stress can be read in Figure 4.8, where 
the row numbers are 917, 918, 919, and 920 and the tensile stresses are 
0.06906 nsi, 0.06905 flsi, 0.06871 llsi, and 0.06872 flsi, while they are not 
exactly uniform but the differences are also smaller than the element in the 
notch region In principle the stress in each element may approach 
uniformity by decreasing the size of the element. This will increase the 
number of the elements and is limited only by the computer memory and 
computation time. The observation of symmetrical model is obvious from 
the stress values of full specimen model and half specimen model. 

B. OBSERVATION OF LINEAR RELATIONSHIP BETWEEN EXTERNAL LOAD 

AND LOCAL TENSILE STRESS 

The tensile stress of the element number 460 (Figure 4.3) of half 
specimen model which subjected to the displacement B.C. value 0.001 is 
0.002475 Msi (Figure 4.9) and the tensile stress of the same element 
which subjected to the displacement B.C. value 0.004 is 0.0112 Msi (Figure 
4.10). The ratio of tensile stress is almost the same as the ratio of 
displacement, the different value bases on the non-uniform of the stress 
around the notch region and can be accepted. When the element number 230 
(Figure 4.7) of the same model, which is far from the notched region the 
ratio of the stresses (Figure 4.11, displacement B.C. value 0.001 and 
Figure 4.12, displacement B.C. value 0.004) and the ratio of the 
displacement are almost the same and the difference is small compare to 
the different value for the element in the notch regioa Therefore, the 
linear relationship between external load and local stress can be accepted 
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whether the element is in notch region or not. Since the element is small 
enough to be satisfy the uniform stress within element, then the uniformity 
and the linearity can be attained. 

C. STATISTICAL SIZE EFFECT RESULTS DISCUSSION 

Equation 2.13 and 2.14, calculated the reliability of the specimen 
model based on statistical size effect in two dimensional, non- 
homogeneous, uniaxial stress. For application in these equation, the 
displacement boundary condition need to be converted to the average stress 
along the boundary. The definition of the average stress in two dimension 
is: 

{P} - (2 ^ ^ '^pq ^ 

where: pq (element number) = 186, 187, 188, 189, 190, 200, 205, 210, 

215, 220, 225, and 230 for half specimen model (Figure 4.7), and 

pq = 163, 164, 165, 166, 167, 168, 171, 174, 177, 180, 183, 186, 
189, 192, 465, 468, 471, 474, 477, 480, 483, and 486 for entire specimen 
(Figure 4.5 and Figure 4.13). 

The values of stress in each element are shown in Table 4.1 and Table 
4.2 respectively. 

The statistical strength of a structure is characterized by the scale 
parameter 3^ in equation 2.13 ( same as BETAA in computer program, 

Appendix A). A change of the value is therefore related to the strength 

change of the structure. The values of 3 e (equation 2.14) for plate with 

notch, where the values of o<, 3. and are 10, 0.1 flsi, and 0.001 
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0.001 respectively are shown in Table 4. 3- A. Table 4. 3-A also shows 
that the variation of the values of which are calculated from different 
values of external load (PBAR) from different values of displacement 
boundary condition (u'). The values of 3^ which calculated from half 

specimen model seem to be higher than 3^ value from the full specimen 

model. The reason is that the number of elements for full specimen is less 
than the number of elements for half specimen model, therefore, the 
accuracy for full specimen model could be lower than that for half 
specimen model. For the purpose of reliability calculation for the 

structure, one can average the 3f values from one set of finite element 
data output from different values of load under the same boundary 
condition, therefore,for instance the average value of 3^ is approximately 

0.030. Table 4. 3-B shows the values of 3^, where the values of <x, and 
are changed to 20, and 0.0001 respectively. The 3^ values are reduced but 
the difference of the average 3^ values is not large, because of the 
compensation of higher o< and lower k^ values. This means that the 

changes of o< and k^ values do not affect much to the 3^ value but o< affects 
the reliability of the structure (equation 2.13). All the reliability 

calculation results for different values of 3g .where the specimen models 
are operated in the same loading range are shown in Appendix B. 
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In the case of plate without notch that shown in Table 4. 3-C and 
Table 4. 3-D, the values of from different values of displacement 
boundary condition are exactly the same. The reason is that the uniformity 
and the linearity affect the accuracy of the calculations. The value of 3^ 

from the plate without notch (0.0439) is higher than from the plate with 
notch (0.028) under the same boundary condition (u' = 0.01), which is a 
reasonable result. 

In addition, the results show that the material scale parameter, 3 is 
converted to 3^ hy the geometry of the structure and the loading boundary 

condition, but it does not depend on magnitude of the load. The results of 
the specimen model without notch confirmed that the uniform stress 
assumption is valid and known result can be recovered. Between results of 
notched specimen model we showed that 3^ remains constant for different 

load magnitude because the uniformity and the linearity conditions are 
attained. As expected we confirmed the notched model is less reliability 
than the plain model under the same loading condition. Given available 
experimental data of a composite material, we can evaluate the scale 
parameter 3 and shape parameter <x. Using these results we demonstrated 
that the procedure of analyzing the effect of a non-uniform state of stress 
obtained from the finite element method can be post processed by 
statistical reliability formulation accounted for size effect, resulting on 
the prediction of the reliability of the structure. 
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Figure 4.1 Element Numbers of Entire Specimen Model at Notch 
Region 
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Figure 4.2 Stress Distribution Values of Entire Specimen Model 



Based on 4-Point Integration, where u' = 0.01", at 
Notch Region 
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Figure 4.4 Stress Distribution Values of Half Specimen Model Based 
on 4- Point Integration, where u' = 0.0r’, at Notch 
Region 
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Figure 4.10 Stress Distribution Values of Half Specimen Model Based 
on 1-Point Integration, where u' = 0.004", at Notch 
Region 
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Figure 4.12 Stress Distribution Values of Half Specimen Model Based 
on l-Point Integration, where u' = 0.004”, at Right 
Boundary Region 
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Figure 4.13 



Element Numbers of Entire Specimen Model at 
Bottom-Right Boundary Region 
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TABLE 4.1 



AVERAGE 5TRE55 FOR HALF SPECIHEN HODEL WITH NOTCH 



ELEMENT 


^ 

<7y , Msi. 


NUMBER 


WIDTH, IN. 


u' = o.oor 


U' = 0.004" 


u' = o.or 


186 


0.057000 


0.009071 


0.036280 


0.090670 


187 


0.057000 


0.009358 


0.037420 


0.090350 


188 


0.057000 


0.009496 


0.037970 


0.094830 


189 


0.057000 


0.009554 


0.038190 


0.095370 


190 


0.057000 


0.009750 


0.038280 


0.095570 


195 


0.024875 


0.005758 


0.035030 


0.087550 


200 


0.024875 


0.008491 


0.033960 


0.084860 


205 


0.024875 


0.008179 


0.032710 


0.081700 


210 


0.024875 


0.007840 


0.031350 


0.078260 


215 


0.024875 


0.007500 


0.029980 


0.074810 


220 


0.024875 


0.007195 


0.028760 


0.071720 


225 


0.024875 


0.006968 


0.027840 


0.069400 


230 


0.024875 


0.006855 


0.027880 


0.068220 




pbar’* 


0.008717 


0.034850 


0.087045 



AveragQ boundary stress 
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TABLE 4. 2 



AVERAGE STRESS FOR FULL SPECIMEN MODEL WITH NOTCH 



ELEMENT 


(?y , Msi. 


NUMBER 


WIDTH, IN. 


u' = o.or 


163, 166 


0.095000 


0.091560 


164, 167 


0.095000 


0.094940 


165, 168 


0.095000 


0.095670 


171, 465 


0.024875 


0.087600 


174, 468 


0.024875 


0.084860 


177, 471 


0.024875 


0.081660 


180, 474 


0.024875 


0.078170 


183, 477 


0.024875 


0.074670 


186, 480 


0.024875 


0.071540 


189, 483 


0.024875 


0.069200 


192, 486 


0.024875 


0.068030 




PBAR** 


0.078025 



Average boundary stress 
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TABLE 4.3-A 

STATISTICAL DATA SUnilARY 
FOR PLATE WITH NOTCH , o< = 10, 3 = 0.1, k^,= 0.001 



MODEL 


HALF 


FULL 




u'=o.oor 


U'=0.004" 


u'=o.or 


U’=0.01" 


pbar’’ 


0.008717 


0.034850 


0.087045 


0.087025 


h 


0.030582 


0.030749 


0.031110 


0.028880 



TABLE 4. 3-B 

STATISTICAL DATA SUMMARY 
FOR PLATE WITH NOTCH, (x = 20. ^ = Oil. k^= 0.0001 



MODEL 


HALF 


FULL 


tisr 


u'=o.oor 


U-0.004" 


u'=o.or 


u'=o.or 


pbar’' 


0.008717 


0.034850 


0.087045 


0.087025 


h 


0.029122 


0.029298 


0.029721 


0.028826 



Average boundary stress 
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TABLE 4.3-C 

STATISTICAL DATA SUMflARY 
FOR PLATE WITHOUT NOTCH. o< = 10, 3 = 0.1, k^= 0.001 



— — WSPUCEMEHT 

UST ' — — 


u’=o.oor 


u'=o.or 


PBAR* 


0.009512 


0.094910 


k 


0.043948 


0.043948 



TABLE 4. 3-D 

STATISTICAL DATA SUriMARY 

FOR PLATE WITHOUT NOTCH . o< = 20, 3 = 0.1, 0.0001 



UST ■ — 


u'=o.oor 


u'=o.or 


pbar'’ 


0.009512 


0.09491 


k 


0.051800 


0.05180 






Average boundary stress 
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APPENDIX A 



STATISTICAL STRUCTURAL RELIABILITY PROGRAH 



******>A****»**«A« *************** »**«***«**********k****A** 
************************************************************ 
PROGRAM: reliability CALCULATIOM OF PuAie MODEL BASED ON 

linearize size effect. NONONIFORM uniaxial stress. 
REHABILITt formulation is expressed in the 
" STANOaPO WEIBOLL FORM 

Data input: output of finite element code name • nikeeo" 

*A**«*«********#ftA**A**R*'*******’R**R*«**4*****A**A****A«A4** 

VARIABLE: KC = MATERIAL CONSTANT FOR PLATE THAT IS 

SUBJECTED TO UNIFORM TENSILE STRESS. IH 
THIS CASE. ASSUME IN EACH ELEMENT THE 
LOCAL TENSILE STRESS IS UNIFORM. 

BETA = SCALE PARAMETER FROM EXPERIMENTAL RESULT 
FROM SEVERAL NUMBERS O'- SAMPLE. 

alpha = shape parameter from expeptmental result 

FROM SEVERAL MUM0ERS OF SAMPLE. 

PBAR : AVERAGE STRESS ALONG BOUOAPY SIDE OF THE 

PLATE THAT subjected TO THE EXTERNAL FORCES. 

IT must be calculated separately after obtaining 
THE STRESS RESULTS FROM "NIKE2D" PROGRAM. 

N£ = number of elements OBTAINING FROM MESH GENERATOR 
PROGRAM NAME "MAZE" WHICH IS THE PROGRAM THAT 
CREATE AN INPUT FILE 70 BE USED IN "NIKEBO'' 

program, also boundary Element numbers can be 
created from “MAZE" PROGRAM. 



DOUBLE PRECISION F . FA .FF . P(. lOOOT . R A( 1 000 ) . BETAA 
OIHENSION XL(2563t .XW<2S6a).SIGYL35b3) 

REAL KC 
byte AL30J 
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ooooooooo 



5 



I 



2 

22 



23 

24 
3 



ryPE*. 'E nter fUE CONTAINING STRESSES DATA INPUT’ 

READC5,5 ) a 
forma T(30AI) 

0 P£N(UNITs 4, F1 LE=A. status* ’OUD’) 

0PENfUNrr=6, FIL£=’SS.0AT’ , STATUS* ’ NEW’ ) 

write ( 6,1) 

F0RMAT(3X, ’DATA OUTPUT FOR L1NEARI7E SUE EFFECT 
*UNIAXIAL TENSIUE STRESS’) 

TYPE*, ’ENTER 1 FOR plate KITH NOTCH 

* OTHERWISE FOR PLATE WITHOUT NOTCH' 

READ»,NC 

IF(NC.EQ.I) THEN 
WRTTE(6,2 ) 
else 

WRITE(6,22) 

END IF 

F0fiMAT(2OX, ’FOR PLATE WITH NOTCH’) 

FORMATCiaT. ’FOR PLATE WITHOUT NOTCH’) 

Type*, ’ENTER the value of DISPLACENENT boundary condition’ 

READ*,U 

WRir£(b,3) u 

type*, ’ENTER I FOR ENTIRE SPECIMEN INPUT 

* otherwise for half specimen INPUT’ 

READ*, NS 

IF(N3.EQ.1 ) then 

WRl IE(6,23) 

ELSE 

wRlTE(b.24) 

END IF ' 

FORHATdSX, ’ENTIRE SPECIMEN MODEL INPUT’) 

F0RMAT( 16X, ’HALF SPECIMEN MODEL INPUT’) 

FORMAT! 7X, 'WHERE THE DISPLACEMENT 8OUN0ART CONDITION IS ’,F8 .U) 
type*, ’ENTER The KC VALUE’ 



PEAD*,KC 

type*, 'ENTER the beta VALUE’ 
READ*, beta 

TYPE*, 'ENTER THE ALPHA VALUE’ 



READ* .ALPHA 
ALPHA2 = ALPHA-f2.0 
FFsO.a 

TYPE*., ’FNTER THE PBAR VALUE’ 
PEAD*,PBAR 

type*, 'ENTER THE NUMSFR op ELEMENTS’ 
P£AD*.NE 



READ element length, WIDTH, ANO' TENSILE STRESS FROM INPUT FILE 



CALCULATE 8 ETAA WHICH IS BETA THAT CHANGED WITH GEOMETRY OF THE 
SPECIMEN MODEL ANO CERTAIN LOADING SOUNOARt CONDITION au f IT IS 
INOEPENOENf OF THE VALUE OF LOAD UNDER TH£ SAME BOUNDARY CONDITION 
DO 200 1-T.NE 

fi£A0(.4,20) XL(.n*'XWtr>.SI&YCI) 

IFISIGYlD.cr.O.O) THEN 
F*(XL(r)*YW(I)/KC**2.0) 

• »(31GT(I)**ALPHA2) 

PF*FF*F 

ELSE 
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non 



FFsfF 
£N0 JF 

200 continue 

IFfNS.EQ.l) THEN 
FFsFF 
EL3£ 

FFsFF*2.0 
END IF 

FA=FF**C-I ./ALPHA2) *P6AR 
BE TAA = BeTA*»( ALPHA /ALPf1A2)-*FA 
WRirE( 6 . 10 ) KC 
10 F0RMAT(8X, 'KC s',Ft3.8) 

WRITEC6.U) beta 
n F0RMAT(6X, ’BETA =• .F13.9) 

WRITE(6.12) alpha 
IP FOR^Af (3x, ’ALPHA =',F13.8) 

WRirE(6,13) PBAR 
13 FOP-'A r (6X, ’PBAR =’,F13.8) 

WRirEl6,lA) NE 
1 a foR'’a r (ax, ■ ne =• ,I4) 

NRIIECb.lS) ALPHA2 
15 F0R«Af(4X, ’ALPHA2 =’,F13.S) 
WRITE(6,25) BETAA 

type*, ■ enter THE FIPSI external LQAD’ 
RE40»,P( 1 ) 

PtMC=P( 1 )/30. 



PELTABUITr CAtCULATIOM BASED ON LINEARIZE SIZE EFFECT. 



1 00 
?0 
2S 
30 



DO 100 J=1,100 

RACJ)=DEXP(,-(P( J1/8ETAA) **ALPHA2) 

Wfttr£(6,30) J,P(J),RA(J) 

PU*1 )=P(J)yPINC 
CONTINUE 

format ( 8X,2F11. Or E« 2.0) 

FoRmaTisX, ’BETAa =*,F13.5) 

FORNArfBX.’T i',IA,3X,'P = ’ ,F 1 3 .8»3X , ’RA =',FI3.3) 

STOP 

END 
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APPENDIX B 



DATA OUTPUT FOR LINEARIZE SIZE EFFECT UNIAXIAL TENSILE STRESS 

data OUTPOT for l-IN£ARrZE SIZE EFFECT UNIAXIAL TENSILE STRESS 
for Plate with notch 

WHERE Th£ displacement BOONOART CONDITION IS O.OUIO 
Half specimen model input 



KC 




a. 


oooioooo 








beta 


- 


0 . 


10040000 








Alpha 




20 . 


oooooooo 








P8AR 


= 


Q. 


00071700 








N£ 


s 


482 










ALPHA2 


r 


22. 


OOOOQOOO 








&ETAA 


= 


0. 


02912212 






J 




1 


P = 


0.00071 700 


RA s 


1 .OOOOOOOO 


J 


- 


2 


P 2 


0 .00000 7S7 


RA = 


1 .OOOOOOOO 


J 


s 


3 


P = 


0.00929013 


RA = 


1 .oooooooo 


J 


s 


'i 


P s 


0.00950'* 70 


Ra 5 


1 .oooooooo 


J 


2 


S 


P 2 


0.00957927 


RA = 


1.00000000 


J 


2 


0 


P 2 


0.01016945 


RA 2 


1 .oooooooo 


J 


2 


7 


P 2 


0.01046040 


RA s 


1 . 0 0 0 0 0 0 0 u 


J 


2 


A 


P 2 


0.')107S097 


RA 2 


1 .oooooooo 


J 


2 


0 


P = 


0.011 )4 153 


RA 2 


1 .oooooooo 


J 


2 


1 0 


P 2 


O.un 3321 4 


RA 2 


1.00000000 


J 


2 


1 1 


P 2 


O.Ol 162267 


RA 2 


1 .oooooooo 


J 


2 


12 


P 2 


O.Ol 191 523 


RA 2 


1 .oooooooo 


J 


2 


1 5 


P = 


0.01220580 


RA 2 


1 .oooooooo 


J 


2 


14 


P 2 


0.01249457 


RA 2 


0 , RRqRRQRR 


J 


2 


IS 


P = 


u. 01278495 


R A 2 


0.99999999 


J 


2 


16 


P 2 


'). 0 1 30 7 5S0 


RA 2 


0 .99999998 


J 


2 


1 7 


P 2 


0.0 1 356607 


RA 2 


0.99999996 


J 


2 


IH 


P 2 


0.01 565665 


RA 2 


0.99999994 


1 


2 


1 9 


0 2 


0.01394720 


RA 2 


0 . 9999999 1 


■J 


2 


20 


P 2 


0.01423777 


2 


0 ,99999085 


J 


2 


21 


P 2 


0.01452835 


RA 2 


0.99999977 


J 


2 


22 


P 2 


0.01481 8<70 


RA 2 


0 ,99999965 


J 


2 


25 


P r 


0.01510947 


PA 2 ' 


0.99999946 


J 


2 


?4 


P s 


0.0 1 5400<»3 


R A 2 


0.9999991 B 


J 


2 


as 


M 2 


0.0 1 5690*30 


R A = 


0.99999877 


J 


2 


26 


P 2 


O.Ol S981 1 7 


RA 2 


0.999998 1 5 


J 


2 


?7 


P 2 


0.01627175 


RA 2 


0 , 999997^5 


J 


2 


h 


P 2 


0.0 1 t^5o23 0 


RA 2 


0 . 9999959 s 


J 


2 


29 


P 2 


0.0 1 685 2^^7 


RA 2 


0 • 99999406 


J 


2 


30 


0 2 


0.01 714545 


RA = 


0.999O913S 


J 


2 


31 


P 2 


0 . 0 1 7 J 5 4 0 0 


R A 2 


0.90998747 


J 


2 


32 


P 2 


O.Ol 772 1S7 


R A 2 


0 ,99990198 


J 


2 


33 


P 2 


0 . 0 1 a 0 1 5 1 5 


RA 2 


0.99907425 


J 


2 


34 


P 2 


0.01850370 


RA 2 


0.99996336 


J 


2 


3S 


P 2 


0.01859627 


RA 2 


0 . 9999481 9 


J 


2 


36 


P 2 


0.01 8386'*3 


RA 2 


0.99992712 


J 


2 


37 


P 2 


0.01917740 


RA 2 


0 .99989804 


J 


2 


3d 


P 2 


0.0l9'ib7o7 


R A 2 


0.999«S806 


J 


2 


39 


P 2 


0.01975855 


RA 2 


0.99980337 


J 


2 


<tO 


P 2 


0.02004910 


RA 2 


0.99972890 


J 


2 


41 


P 2 


0 ,020 5.5967 


RA 2 


0.99962797 


J 


s 


a2 


P 2 


0.0206 5 )?5 


RA r 


0.999491 75 


J 


2 


45 


P 2 


O.OPO 4^-)80 


PA 2 


0.99950870 


J 


2 


U4 


P 2 


0.02) 21 1 57 


PA 2 


0.9990657 5 


J 


= 


4S 


P 2 


0 . 0 2 1 5 0 1 9 5 


R A 2 


0.99075722 


J 


2 


46 


L> - 


0. TP 1 79 25 ) 


PA 2 


0,99830578 


J 


= 


47 


P 2 


0. 022^^507 


RA 2 


O. 99773 O 6 I 


J 


2 


48 


? 2 


'^.O22 5 7 5o5 


RA 2 


0.99697558 


J 


2 


49 


0 - 


0 .02266420 


RA 2 


0.99898076 
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Data output for LINeARlZe SHE EFFeCT uniaxial tensile stress 
FOR PLATE VXITH NOTCH 

WHERE THE OISPLACEMCNT 90UM0ARY CONDITION 18 0.0100 

half specimen model input 



KC 


- 


0, 


0001 0000 








BETA 


z 


0. 


10000000 








Alpha 


z 


20. 


00000000 








P8AR 


= 


0. 


08704300 








M£ 


z 


482 










ALPHAS 


z 


22. 


00000000 








8£TaA 


z 


0. 


02972093 








J 


z 


1 


P s 


0. J0871 700 


PA z 


1 .00000000 


J 


z 


2 


P r 


0.;)09n07S7 


PA Z 


1 .00000000 


J 


z 


s 


P = 


0 . 00<^Po^ 1 z 


RA z 


1 .00000000 


J 


s 


a 


P = 


0.009S8970 


RA z 


1 .00000000 


J 


z 


5 


P s 


0.00987927 


RA z 


1 .00000000 


J 


s 


o 


P r 


0.0l0lo983 


RA r 


1 .00000000 


J 


z 


7 


P = 


o.oioa^oao 


RA z 


1. 00000 0 00 


J 


z 


9 


P = 


0.01075)07 


PA z 


1 .00000000 


J 


z 


9 


P r 


0,01 1 OiitS5 


RA z 


1 .00000000 


J 


z 


10 


P - 


0 . 0 1 1 3 3 2 1 0 


RA r 


1 .00000000 


J 


z 


1 1 


P z 


0,011 822o7 


PA z 


1 .00000000 


J 


z 


12 


P z 


O.Ol 191523 


RA r 


1 .00000000 


J 


z 


15 


P z 


0.01220350 


RA z 


1,00000000 


J 


z 


1 a 


P z 


0.01209457 


RA z 


0.99999999 


J 


z 


IS 


o z 


0.012M495 


RA r 


0,69966696 


J 


z 


lb 


P z 


0.01 307550 


PA r 


0.99996999 


J 


z 


17 


p = 


0.01 536507 


RA z 


0.99999998 


J 


z 


la 


P z 


0.01 3SSbb5 


RA z 


0.99999996 




z 


19 


p = 


0.01 39 '4 72“ 


PA r 


f) ,66969094 


J 


z 


20 


p = 


0.01423777 


PA z 


0,9999669 1 


J 


z 


21 


o z 


0.01052‘^55 


RA r 


0,99999985 


J 


z 


2?. 


p z 


0.01081890 


RA z 


0,9^999978 


J 


z 


23 


P z 


0.015109/17 


RA z 


0 ,999996(36 


J 


z 


2^ 


P z 


0.01500003 


R A z 


0 , 9660QQa« 


J 


z 


25 


P s 


O.Ol5o9 jbO 


RA z 


0.99999921 


.1 


z 


26 


3 z 


0.01598117 


RA z 


0.99999882 


J 


z 


27 


P z 


0.016271 75 


PA z 


0 , 06669^25 


J 


z 


2S 


P z 


0.0l65o250 


RA z 


0 ,99999741 


J 


z 


29 


P z 


O.Ul6«5287 


RA r 


0 ,99999620 


J 


z 


30 


o z 


0.01 714303 


PA z 


0 ,9096944 7 


J 


z 


31 


P z 


0.01745/400 


PA z 


0,99999200 


J 


z 


3? 


P z 


0.01 772457 


RA z 


0.99990846 


J 


z 


33 


P z 


0.01801513 


PA r 


0.99998353 


J 


z 


3^ 


P z 


0.0 1 830570 


PA z 


0.99967658 


J 


z 


35 


P z 


0.01^59627 


PA z 


0.99966689 


J 


z 


3b 


P z 


0.018^8683 


R A z 


0.9996S343 


J 


z 


37 


P z 


0.01917740 


PA z 


0.69993484 


J 


z 


38 


P z 


0.0194o797 


PA z 


0.9999092^^ 


J 


z 


39 


P z 


0.01975353 


RA z 


0.999«7434 


J 


z 


UO 


P z 


0.02004910 


RA z 


0.99982675 


J 


z 


u 1 


p z 


0.02033967 


PA z 


0.99976224 


J 


z 


ae 


P z 


0.02063023 


RA z 


0.9666751^ 


J 


z 


az 


P z 


0.02092980 


PA z 


0.99955817 


J 


z 


uu 


P z 


0.02121 1 37 


fVA z 


0 .96640 1 58 


J 


z 


as 


P z 


0.02150193 


PA z 


0.96019284 


J 


z 


Ub 


P z 


O.o2 1 79 250 


PA z 


0,96891570 


J 


z 


(i7 


P z 


0.oP20^S07 


RA z 


0,99854916 


J 


z 


ae 


P Z 


0.0P2373o3 


PA z 


0.96806620 


.1 


z 


49 


P r 


0.0226o42'‘ 


R A z 


0.09743222 



65 



data output F0« linearize size EPFECr uniaxial tensile stress 

FOR PLATE WITH NOTCH 

WHERE THE DISPLACEMENT BOUMDARt COnOITION IS 0.00<»0 

Half specimen model input 



KC 


r 


0. 


OOOiOOOO 








beta 


- 


0. 


10000000 








ALPHA 


- 


20. 


oooooooo 








P6AR 


= 


0. 


03485000 








N£ 


z 


462 










ALPHA2 


= 


22. 


OOOOOOOO 








betaa 


= 


0. 


02929771 








J 


z 


1 


P = 


). 008 71700 


PA 3 


1 .oooooooo 


J 


z 


2 


P z 


0. 00900757 


PA 3 


1 . 0 0 0 0 0 0 0 0 


J 


z 


3 


P z 


0. J0OP931 3 


PA 3 


1 .oooooooo 


J 


z 


a 


P z 


0.0095 ^"^70 


PA 3 


1 .oooooooo 


J 


: 


5 


P 3 


0.009M7927 


PA 3 


1 .oooooooo 


J 


z 


6 


P Z 


O.ul 016963 


PA 3 


1 .oooooooo 


J 


z 


7 


p z 


0 . 0 1 0 4 6 0 4 0 


PA 3 


1 .oooooooo 


J 


z 




P z 


0.01075097 


PA 3 


1 .oooooooo 


J 


z 


9 


P z 


0.01104155 


PA 3 


1 .oooooooo 


J 


z 


10 


P 3 


0.01153210 


RA 3 


1 .OOOOOOOO 


J 


z 


1 1 


P 3 


0.1)1 1 6226/ 


PA 3 


1.00000000 


J 


z 


12 


P 3 


0.01191 523 


PA 3 


1 .oooooooo 


J 


z 


1 3 


P Z 


0.012?03«'» 


RA 3 


1 .oooooooo 


J 


z 


1 4- 


P 3 


0.v)1249*j57 


PA 3 


0.99999999 


J 


z 


15 


P 3 


0 .0 1 27H49 3 


PA 3 


0 . pqpQPQPP 


J 


z 


1 6 


P - 


O.Ol 307550 


PA 3 


0.99999998 


J 


z 


1 7 


P 3 


O.Ol 336607 


PA 3 


0.99999997 


J 


z 


18 


P Z 


0.01 365663 


PA 3 


0 .9994999*^ 


J 


z 


1 9 


P 3 


0.01 594720 


PA 3 


0 .63999002 


J 


z 


20 


P z 


0 . 0 1 a ? 5 7 7 7 


OA 3 


0.99990987 


J 


z 


21 


P z 


0.01452833 


PA 3 


0.99999980 


J 


z 




P z 


. p 1 a 8 1 


PA 3 


0,09009969 


J 


z 


2 S 


P 3 


0.01510947 


PA 3 


0 .09090053 


J 


z 


2'J 


P z 


0 . 0 1 5 1 0 0 0 3 


Pa 3 


0 .99090P2« 


I 


z 




P 3 


Q .0 1 56^004 


f> 

II 


0 .99990892 


J 


z 


?h 


P z 


0.fil59Hi 1 / 


PA 3 


0.99990836 


J 


z 


?7 


P 3 


).0lr>271 73 


PA 3 


() ,99000759 


J 


z 


2« 


P 3 


0.0lh5-*230 


PA 3 


0 .99000645 


J 


z 


5 0 


P z 


0 . J ! 6852"* 7 


RA 3 


0,09000479 


J 


Z 


50 


P 3 


).0l 710 54 3 


PA 3 


0 ,99999242 


J 


z 


31 


P Z 


0 . 0 1 7 4 3 4 0 0 


PA 3 


0.09998903 


J 


z 


32 


P 3 


O.Ol 772457 


RA 3 


0.9999H421 


J 


z 


33 


P 3 


0.01801513 


PA 3 


0.99097742 


J 


z 


5a 


P Z 


0.01 A30570 


PA 3 


0.99996790 


J 


z 


35 


P 3 


0.01859627 


PA 3 


0 . O9OO«5/i50 


J 


z 


5fe 


P 3 


0 . j 1 6 8 6 8 3 


PA 3 


0.99993615 


J 


z 


37 


P 3 


0.0 191 77 40 


II 

< 

Q 


0.99991067 


J 


z 


58 


P 3 


0.0104^797 


PA 3 


0.99067564 


J 


z 


TP 


P 3 


0.01975853 


RA 3 


0.99982772 


J 


z 


uo 


P 3 


0.02004910 


PA 3 


0.99076248 


J 


z 


'4 1 


P Z 


0.02053967 


□A 3 


0.99967405 


J 


z 


42 


P Z 


0.u2ud3)23 


PA 3 


0.99955470 


J 


z 


4 3 


P z 


0. 02092060 


PA 3 


0.99039431 


J 


z 


4 4 


P 3 


0 . 0 2 1 2 1 1 3 7 


PA 3 


0.9991 7967 


J 


z 


45 


P 3 


0.02150195 


PA 3 


0.99689357 


J 


z 


4^ 


P Z 


0.021 79250 


PA Z 


0.99851 375 


J 


z 


07 


P 3 


0.O22085O7 


RA 3 


0.9980 1147 


J 


z 


aa 


P Z 


0 . )22 37 3o3 


PA 3 


0.99734976 


J 


z 


40 


P 3 


0 .oP2664P<* 


PA 3 


0.996481 50 



66 



qata output for lIN£ariz£ Size Efpect uniaxial tcnsilp stress 
for plate viithoht notch 

WHERE THE displacement SOONOART CONDITION IS 0.0010 
ENTIRE specimen MOOEL INPUT 



Kc 


: 


0. 


0001 0000 








BETA 


= 


0, 


lOOOOOOO 








ALPHA 




ao. 


OOOOOOOO 








PBAR 


= 


0. 


00951200 








NE 


s 


200 










ALPHA2 


z 


22. 


OOOOOOOO 








8ETAA 


z 


0. 


O51B0000 








J 


z 


1 


p = 


O.OOH71 700 


z 


1 .OOOOOOOO 


J 


z 


2 


P r 


0.00900757 


PA s 


1 .00000000 


J 


z 


3 


P s 


O.OOOPOAi S 


RA z 


1 .OOOOOOOO 


J 


z 


a 


P = 


0.0095»*^7<> 


RA r 


1 .OOOOOOOO 


J 


z 


•5 


P = 


0.O09M7927 


RA r 


1 .OOOOOOOO 


J 


z 


b 


P 5 


0,0101g9«^ 


RA r 


1 .OOOOOOOO 


J 


z 


7 


P = 


0.0 1 UUbOaO 


RA z 


1 .OOOOOOOO 


J 


z 


A 


f’ 5 


0.0 1 075097 


RA s 


1 .OOOOOOOO 


J 


z 


9 


P = 


0.01 ioai53 


RA z 


1 .1)000 0000 


J 


z 


1 0 


P = 


0 . 0 1 1 3 3 2 1 0 


RA r 


1 .OOOOOOOO 


J 


z 


1 1 


P = 


0 .0 1 1 b2267 


RA r 


1 .OOOOOOOO 


J 


z 


12 


P = 


0.01 1 91 523 


RA r 


1 .OOOOOOOO 


J 


z 


13 


3 = 


0.01220580 


RA r 


1 .OOOOOOOO 


J 


z 


ly 


P r 


0.0l2'49a37 


PA z 


1 .OOOOOOOO 


J 


z 


15 


P = 


0.01278995 


RA z 


1 .OOOOOOOO 


J 


z 


lb 


P = 


0.0 1 5 0 7 550 


RA z 


1 .OOOOOOOO 


J 


z 


1 7 


P = 


0.Ol53bb07 


RA z 


1 .OOOOOOOO 


J 


z 


Id 


P = 


0.015b5bb^ 


RA z 


1 .OOOOOOOO 


J 


z 




P = 


0.01 30a720 


RA z 


1 , 0 0 0 0 0 0 0 0 


J 


z 


20 


P = 


0.01023777 


RA z 


1 .OOOOOOOO 


J 


z 


21 


P z 


O.CM 45285 5 


RA r 


1 .OOOOOOOO 


. J 


z 


bP 


P = 


0 . ) 1 4 3 1 HO n 


RA z 


1 .ooooonoo 


1 


z 


2 ^ 


P • 


i).iM5l099 7 


R A z 


1 .OOOOOwOO 


J 


z 


2^^ 


P = 


0.01 540005 


RA z 


1 . 0 0 0 0 0 0 0 0 


J 


z 


25 


P = 


O.'USbPObt- 


R A z 


1 . 0 0 0 4 0 :) U 0 


J 


z 


2b 


o - 


0.. 115981 1 7 


R A z 


1 .OOOOOOOO 


J 


z 


27 


O ; 


0 . 0 1 o2 7 1 7 3 


P A z 


1 .OOOOOOOO 


J 


z 


26 


p = 


u . 0 1 6 5 o 2 5 0 


RA = 


1 .ooooooon 


J 


z 


29 


p = 


0 . >1 1 o h5 25 7 


PA z 


r.oooooooo 


J 


z 


30 


P r 


0. M 7 14 30 5 


R A = 


1 .OOOOOOOO 


J 


z 


31 


p = 


0 . 0 1 7 4 5 u 0 0 


PA z 


1 . 0 n 0 u 0 0 0 0 


J 


z 


32 


P T 


O.ol 772457 


•^A z 


1 . 0 0 0 0 0 0 0 0 


J 


z 


35 


P r 


0.0 1 HO 1 51 3 


R A z 


1 .OOOOOOOO 


J 


z 


3^ 


P s 


0.01030570 


RA z 


1.00000000 


J 


z 


3‘5 


P = 


0.0tM59bP7 


RA z 


1 .OOOOOOOO 


J 


z 


3rt 


P = 


■). )1H88S55 


R A z 


1 .000000 00 


J 


z 


37 


p s 


0.01 PI 7 740 


R A z 


1 .OOOOOOOO 


J 


z 


38 


p = 


0. 01^4^797 


PA z 


1 .OOO'IOOOO 


J 


z 


39 


p = 


0.0 19 75H5.3 


RA z 


1 .OOOOOOOO 


J 


z 


^0 


P s 


0 . )0OOa9 1 0 


RA z 


1 .OOOOOOOO 


J 


z 


U.) 


p r 


0.02033967 


R A z 


1.00000000 


J 


z 


y-2 


P s 


0.020‘,5025 


PA z 


1 .OOOOOOOO 


J 


z 


a3 


o r 


0.0209?;)HO 


R A z 


1.090 00000 


J 


z 




P = 


0 . 0 a 1 a 1 1 5 7 


R A z 


1 .OOOOOOOO 


J 


z 




P = 


o.oai 501 ^5 


PA z 


1 .OOOOOOOO 


J 


z 


a6 


P = 


0.021 792h0 


RA z 


o'. 99999P99 


J 


z 


(il 


P = 


0 . 0 P P 0 H 5 0 7 


RA z 


0 . :?4999P9P 


J 


z 


aS 


P = 


0. 022 37 363 


PA z 


0 , pqogqogo 


J 


z 


9.9 


P = 


0.022*'^ I2f' 


R A z 


0 . 999RQ99P 
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data output for linearize size effect uniaxial tensile stress 

FOR PLATE WITHOUT NOTCH 

WHERE THE DISPLACEMENT BOUNDARY CONDITION IS 0.0100 
ENTIRE SPECIMEN MODEL INPUT 



KC 


z 


0. 


0001 0000 








BETA 


z 


0, 


tooooooo 








ALPHA 


z 


20. 


oooooooo 








PBAB 


z 


0. 


09AB1000 








Ht 


z 


200 










ALPHA2 


z 


22. 


oooooooo 








0£TAA 


s 


0. 


05180009 








J 


z 


I 


p = 


9.00b7l 7<H) 


PA r 


1 .OOOOOOOO 


J 


z 


2 


p = 


0.009ri07S7 


z 


1 . 0 0 0 0 0 0 f) 0 


J 


z 


5 


p = 


0. )0P29o I 5 


PA z 


1 , 0 0 0 0 0 0 0 0 


J 


z 




P z 


0.')O9oR'^7n 


PA z 


1 .oooooooo 


J 


z 


5 


P s 


0,009S7P^ 7 


RA z 


1 .oooooooo 


J 


z 


6 


P z 


0.01016985 


RA z 


1 .oooooooo 


J 


z 


7 


P z 


0.0 1046040 


RA z 


1 .oooooooo 


J 


z 




p z 


0.0 10 7509 7 


RA z 


1 .oooooooo 


J 


z 


9 


P z 


0 . t> n :) 4 IS 5 


PA z 


1 . 0 0 0 0 0 0 0 0 


J 


z 




P z 


0. 01 1 55210 


RA z 


1 .oooooooo 


J 


z 


I I 


P z 


O.Ol 162267 


PA z 


1 .oooooooo 


J 


z 


I2 


P r 


0.01191525 


PA z 


1 .oooooooo 


J 


z 


15 


P s 


0.012205^0 


RA z 


1. 00 0 00000 


J 


z 


I'l 


P z 


0.01 2404 57 


RA z 


1 .oooooooo 


J 


z 


IS 


P z 


0.0 1^78403 


PA r 


1 .oooooooo 


J 


z 


lb 


P z 


0.01507530 


PA z 


1 .oooooooo 


J 


z 


1 7 


p z 


0.0 1 55ob07 


PA z 


1 .OOOOOOUO 


J 


z 


IS 


P z 


0. J 1 56S665 


PA z 


1 . oooooooo 


.l‘ 


z 




o - 


0.01 594720 


PA z 


1 .oooooooo 


J 


z 


?o 


P 2 


0 . 0 1 4 a 5 7 7 7 


PA z 


1 . 0 0 0 0 0 0 0 0 


J 


z 


? 1 


P z 


0.;)1 4 52^5^ 


PA z 


1 .oooooooo 


J 


z 


?2 


P r 


0.01451 -^90 


PA Z 


1 .oooooooo 


J 


z 


?5 


P z 


0 . iM S 1 0 4a 7 


PA z 


1 .oooooooo 


J 


z 


2'4 


P z 


0 • 0 1 S 4 0 4 0 5 


PA z 


1,00000000 


J 


z 




P z 


0,0 1 S690*,n 


PA z 


1 .oooooooo 


J 


z 




P z 


0 . 4 t S9P 1 1 7 


PA z 


1 .oooooooo 


J. 


z 


a? 


p z 


O..M&2 7 1 75 


PA r 


1 .O-IOOoOOO 


J 


z 


28 


P z 


0 . ' M 5 6 2 5 *. 


RA z 


1 .oooooooo 


J 


z 


PQ 


P z 


0,'Mft=^S247 


RA z 


1 , oooooooo 


J 


z 


3 0 


P r 


n . 0 1.7 1 '4 5 15 


PA z 


1 .oooooooo 


J 


z 


51 


p z 


0 , 1 7 4 54 )4 


PA z 


1 .OOoOOOOO 


J 


z 


32 


P z 


0.01 772 iS7 


PA z 


! .OOOOuOOl) 


J 


z 


33 


P z 


0.01801515 


PA r 


1 .OOOOOOOO 


J 


z 


3a 


p z 


0.01H50570 


PA z 


1 .oooooooo 


J 


z 


iS 


P z 


0.01 


P A z 


1 .oooooooo 


J 


z 


3b 


P z 


0.01 888S«5 


PA z 


1 . OOOOMOOO 


J 


z 


57 


P Z 


0.01 Qi 7 7-40 


PA r 


1 .oooooooo 


J 


z 




P z 


0 , 1 o M ^ 7 0 7 


PA z 


1 . 0 0 0 0 0 0 0 0 


J 


z 




p z 


0.0 1073^^3 5 


PA z 


1 , oooooooo 


J 


z 


uo 


c z 


0.02004410 


? A z 


1 , 0 0 0 0 0 0 0 0 


J 


z 


ai 


P Z 


0,0a055967 


PA z 


1 .0001)0000 


J 


z 




o z 


0.4a065')25 


PA z 


1 . oooooooo 


J 


z 


45 


P r 


0.0209^ >a.) 


PA z 


1 . oOOOOOOO 


J 


z 


ua 


p z 


0 . 0 9 1 a 1 1 5 7 


RO z 


1 ,000-)000:) 


J 


z 


45 


P z 


0,i»2lS01 >5. 


RA z 


1 . OOOOOOUO 


J 


z 


4^ 


P r 


O.OPl 799S'>’ 


PA z 


0 , Q:)i?f)9Qgo 


J 


z 


a 7 


P z 


J.'»aa085o7 


PA z 


0 , qqgzg^QO 


.J 


r 


48 


P z 


0 ,092 5 7 50 5 


PA z 


0 , Q44P96Q9 


} 


: 


a9 


P z 


0.aa2o^4Pn 


PA z 


0 , 9P9OQ9P0 
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DATA OOrpoT FOR LlNEARUE SIZE EFFECT UNIAXIAL TENSILE STRESS 
FOR PLATE WITH NOTCH 

WHERE the D I SPUAC EWENT BOUNDAR'f CONOITIOH IS 0.0100 
ENTIRE SPECIMEN MODEL INPUT 



KC 


z 


0. 


OOOlOQOO 








BETA 


= 


0« 


10000QOO 








ALPHA 


z 


20, 


00000000 








P8AR 


Z 


0. 


08702500 








NE 


z 


392 










ALPMA2 


z 


22. 


OOOOOOOO 








HETAA 


Z 


0. 


00288262 






J 


z 


1 


P z 


0.00P71 700 


PA z 


1 .0 0000000 


J 


z 


2 


P = 


0,0090 07S7 


RA z 


1 . 0 0 0 0 0 0 0 n 


J 


z 


3 


P Z 


0,0092951 3 


RA z 


1 .oooooooo 


J 


z 


4 


p = 


0,00958570 


RA z 


i .oooooooo 


J 


z 


5 


P z 


0,OOPH7927 ' 


RA z 


1 .oooooooo 


J 


z 


6 


P z 


0,01016983 


RA z 


1 .oooooooo 


J 


z 


7 


P Z 


0.0104b0'40 


RA z 


1 .oooooooo 


J 


z 


e 


P z 


0.0 1 0 75097 


RA z 


1 .oooooooo 


J 


z 


9 


P Z 


0.01104153 


RA z 


1 .oooooooo 


J 


z 


10 


P z 


0,011 3321 0 


RA z 


1 .oooooooo 


J 


z 


1 1 


P Z 


0,01 162267 


RA z 


1 .oooooooo 


} 


z 


12 


P Z 


0,01 1 91 523 


RA z 


1 .oooooooo 


J 


z 


1 3 


D Z 


0,01220580 


RA z 


0.99999999. 


J 


z 


10 


P Z 


0. 01349^57 


PA r 


0 ,99909990 


J 


z 


1 5 


0 z 


0.01276 19^ 


RA z 


0.99999998 


J 


z 


16 


P z 


0,01 307550 


RA z 


0.99999997 


J 


z 


1 7 


P Z 


0.01536607 


R A z 


0.99999095 


J 


z 


IB 


? Z 


0.01 365663 


P z 


0.99909992 


J 


z 


t ^ 


0 z 


n,in 594720 


R A z 


0 .99999988 


J 


z 


20 


3 z 


0.01423777 


RA z 


0.99999981 


J 


z 


21 


P Z 


0.01452835 


RA z 


0,99999071 


J 


z ' 


22 


D - 


0 . 01 1 590 


R A z 


0.9990995/1 


J 


z 


25 


P z 


0.0 15 to 9o 7 


PA z 


0,99999930 


J 


z 




P Z 


0,0 1 5 40003 


R A z 


u . 99999894 


J 


z 


25 


c - 


0 . 0 1 569 0hfi 


R A Z 


0 ^ 99QO98a0 


J 


z 


26 


0 z 


0.015961 1 7 


R A z 


0.99999760 


J 


z 


27 


P z 


0.016271 73 


RA z 


0,99999645 


J 


z 




P z 


0.0 1 656230 


RA z 


Q , 939qq/| 7u 


J 


z 


2^ 


P Z 


0.(il6H5.2*"7 


PA z 


0.99990228 


J 


z 


iO 


P Z 


0.0 1 71 4 sa 3 


RA z 


0,99998876 


J 


z 


31 


D - 


0.01 743400 


RA z 


0.99998573 


J 


z 


32 


P z 


0,01772457 


R A r 


0,99997660 


J 


z 


33 


P z 


0 . 01 H 0 1 5 1 3 


RA z 


0.99996653 


J 


z 


50 


P z 


0.0 18 S05 70 


RA z 


0.93905241 


J 


z 


35 


P z 


n.01'859o27 


RA z 


0,09993271 


J 


z 


3b 


P z 


0.018«8683 


R A z 


0.99900536 


J 


z 


37 


P Z 


0.0191 7740 


II 

<i 


0.99988759 


J 


z 




P z 


0,01 94o797 


RA z 


0.9998 1 567 


J 


z 


30 


P Z 


0.01 9 75*153 


PA z 


0,99974465 


J 


z 


un 


D - 


0,0200491 w 


PA z 


0.99064795 


I 


z 


at 


P z 


0.020 3396 7 


R A z 


0.99951688 


J 


z 


a2 


P z 


0.02063023 


PA z 


0.99934001 


J 


z 


a3 


P z 


0.02092080 


PA z 


0.99910233 


J 


z 


ua 


P z 


0.02121157 


PA z 


0,99878427 


J 


z 


as 


P z 


0.02150193 


PA z 


0,998560 39 


J 


z 


a6 


P z 


0.021 79250 


PA z 


0.99779774 


J 


z 


a7 


□ - 


0, )2?f’^507 


K A Z 


4.99705385 


J 


z 


48 


p - 


0.oP23736^ 


PA z 


0.996074 08 


1 


z 


49 


z 


0 . .12266 xP ‘ 


- A z 


0 . 0 9 .'1 7 8 8 7 0 


I 


z 


50 


- z 


.m2P9'5 i7 / 


; z 


0,O93l08Ml 
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0 *r» OUTPUT FOR LTMEARtZe SIZE EPPECT UNIAXIAL T£NSIl£ STRESS 
FOR PLATE WITH NOTCH 

WHERE THE OrSPLACEMcNT SOONOABY CONDITION IS 0.0010 
half SPEClMeM NOOEL INPUT 



KC 


= 


0. 


ootooooo 










9^r^ 


s 


0. 


10000000 












: 


to. 


oooooooo 










PBAR 


z 


0. 


00871700 










fiE 


= 


482 












ACPHA3 


z 


12. 


OOOOOOOO 










5ZTAA 


z 


0. 


03050189 










J 


z 


1 


P = 


0.')<^67l 700 


RA 


z 


0 .R9qqqq7 1 


J 


z 


2 


o s 


0.00000757 


RA 


z 


0.qqq9qQ57 


J 


z 


5 


P = 


o.ooopq^l 3 


RA 


z 


0.99999938 


J 


z 


u 


P = 


0 .uoqSA5 70 


RA 


z 


0 . 99999010 


J 


z 


5 


P = 


0.0nq97P27 


PA 


z 


0.99999671 


J 


z 


6 


P = 


O.Ol t) 16PH3 


RA 


z 


0.9999981 7 


J 


z 


7 


P = 


0 . 0 1 0 0 b 0 a 0 


RA 


z 


0.99999 744 


J 


z 


A 


P X 


0.01m7«5097 


®A 


z 


0 . 99999ba4 


J 


z 


q 


P = 


0 . 0 1 1 0 q 1 S 5 


RA 


z 


0.99999509 


J 


z 


10 


o Z 


0.01153210 


RA 


z 


0.99999350 


J 


z 


1 1 


p = 


O.Ol 1 b2267 


RA 


z 


0.99999092 


J 


z 


12 


p = 


0.01 1 PI 525 


RA 


z 


0.99998779 


J 


z 


1 5 


o = 


0.01220390 


RA 


r 


0.99998369 


J 


z 




P s 


0.01 ?'iqM37 


RA 


z 


0.99997637 


J 


z 


15 


P 2 


0.O1278UQ3 


RA 


z 


0.99997150 


J 


z 


1 h 


P = 


0.0 1 50 7550 


PA 


z 


0.99996268 


J 


z 


17 


P = 


0 .n 1 33ot)07 


RA 


r 


0.99995142 


J 


z 


l« 


P = 


0.01 3b59o3 


RA 


z 


0.99993712 


J 


z 




P = 


•) , 0 1 57'J 7PO 


RA 


z 


0.99991 904 


J 


z 


20 


P = 


0.0) U23777 


PA 


z 


0.90989632 


J 


z 


21 


P = 


0.0PI52‘<33 


PA 


z 


0.99986787 


J 


z 


22 


P = 


0 .0 1 4^^ 1 


R A 


z 


0 .999 8 524 '4 


J 


z 


23 


P z 


0.01510007 


RA 


z 


0.99978847 


J 


z 


2'J 


p = 


0 . 0 1 5 '4 0 0 0 3 


RA 


z 


0.99975415 


J 


z 


25 


p = 


0.01 SbOObP 


RA 


z 


0.99O5t,752 


J 


z 


2^ 


p = 


0.'M5PPl W 


P A 


z 


0.99958539 


J 


z 


? 7 


p = 


0.0ln?7l 73 


RA 


z 


0.99948554 


J 


z 


2A 


o X 


0.0 1 b5b23^ 


PA 


z 


0.9993(5359 


J 


z 




p = 


0.01^5^5297 


PA 


z 


0.99921595 


J 


z 


30 


p = 


0.01714503 


RA 


z 


0.99903752 


J 


z 


31 


p = 


0 . 0 1 7 4 3 4 0 0 


RA 


z 


0.99882257 


J 


z 


32 


p = 


0.01 772'4S7 


RA 


z 


0.R9856445 


J 


z 


33 


P z 


0.01801515 


RA 


z 


0.99825539 


J 


z 


3M 


p = 


0.01830570 


PA 


z 


0.99786649 


J 


z 


35 


P r 


0.0 185P927 


RA 


z 


0.99744740 


J 


z 


Sh 


P z 


0.01^^589*^3 


RA 


z 


0. 99692624 


J 


z 


37 


P z 


O.OlPl 7 7 0 0 


RA 


z 


0.99630956 


J 


z 


38 


P z 


O.OIP lb707 


RA 


z 


0.995581 10 


J 


z 


55 


P z 


0.01975^53 


PA 


z 


0.99472362 


J 


z 


ao 


o z 


0 . 02O0U0 1 0 


RA 


z 


0.99371655 


J 


z 


<4 l 


p = 


O.OPO 33P87 


PA 


z 


0.99253675 


J 


z 


a2 


o - 


0.0?0b3023 


RA 


z 


0.991 15800 


J 


z 


U5 


P Z 


0 . 0 2 0 0 2 0 8 0 


RA 


z 


0.98955067 


J 


z 


aa 


■p = 


0.0?) 2+1 57 


PA 


z 


0.-9S7681 35 


J 


z 


q.5 


P z 


0.021501P3 


PA 


z 


0.98551251 


J 


z 


at» 


p = 


O.OPl 70250 


RA 


z 


0.98300210 


J 


z 


(17 


p z 


0.02208507 


RA 


z 


0.98010318 


J 


z 


(Ih 


P z 


0.022 573o5 


RA 


z 


0.97676549 


U 


z 


aq 


r • 


0.022bn 420 


V A 


z 


0.97292518 
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data output for linearize size effect uniaxial tensile stress 

FOR PLATE WITH NOTCH 

Where the displacement boundary condition is o.ooao 



HALF specimen WOEL InPoT 



KC 


= 


0. 


oolooooo 








BETA 


= 


0. 


iOOOOOOO 








alpha 




to. 


OOOQOOOO 








PBAR 


= 


0. 


03465000 






ME 


r 


462 










Alphae 


= 


12. 


00000000 








BETAA 


= 


0, 


05074947 








J 


s 


1 


P s 


0.0')e71 700 


RA z 


0.99999973 


J 


z 


2 


P = 


0.0()9'J0757 


94 z 


0. 99999960 


J 


z 


5 


P = 


1.009294 1 5 


RA z 


0.99999902 


J 


z 


a 


P s 


0.0095H4 70 


PA z 


^qg:?99q \ $ 


J 


z 


5 


p = 


0.009*^7927 


RA r 


0.99999879 


J 


z 


b 


P s 


0.Ul01b983 


RA z 


0.99999829 


J 


z 


7 


P S 


O.OlOOoOaO 


RA r 


0.99999760 


J 


z 




o z 


0.01O75097 


RA z 


0.q9999666 


J 


z 


9 


P z 


0.011 04153 


RA z 


0,99999540 


J 


z 


t 0 


P = 


O.Ol 1 33210 


RA z 


0 , qoqoq ^72 


J 


z 


1 1 


P z 


0.011 622b7 


RA z 


0.99999150 


J 


z 


1? 


P s 


0.01191323 


RA z 


0.99998856 


J 


z 


1 i 


P S 


0.01220 580 


PA z 


0 . 99993/473 


J 


z 


1 4 


P s 


0.012494S7 


PA z 


0.99997975 


J 


z 


IS 


P r 


0.01 2 73*193 


RA z 


0,99997331 


J 


z 


1 6 


P = 


0. 01 307550 


RA z 


0.99996505 


J 


z. 


17 


P = 


0.01 336607 


RA r 


0.99995450 


J 


z 


18 


P z 


0.01355663 


RA z 


0.999941 1 1 


J 


z 




P z 


0.01 394720 


PA z 


0.99992418 


J 


z 


20 


P z 


9.01423777 


PA Z 


0 ,99900290 


J 


z 


21 


p z 


0.01452533 


RA r 


0.99987626 


J 


z 


22 


p = 


o.oia^tspn 


RA z 


0.99984307 


J 


z 


23 


p = 


0.01510947 


RA z 


0 .999.8 0 1 90 


J 


z 


24 


P z 


0 . 0 1 5 'J 0 0 f» 3 


R A r 


0.99975 1 <1 3 


J 


z 


25 


P z 


0,015590o0 


PA r 


0.99968843 


J 


z 


5r» 


P z 


0.0 ISP31 1 7 


PA z 


0.99961 1 70 


J 


z 


27 


? = 


0 . 1 b 2 7 1 7 3 


RA z 


0.9.9651800 


J 


z 


2a 


P z 


0.01o5o23) 


R A z 


0.9994 0397 


J 


z 


29 


P z 


O.Olb'^GPd/ 


R A z 


0.99926570 


J 


z 


30 


p z 


0 . l) 1 7 1 4 3 4 3 


R A z 


0.9990985^ 


J 


z 


31 


p = 


0.01 743400 


RA z 


0.99389726 


J 


z 


32 


P z 


O.ni 7 7205 7 


R A z 


0.99865549 


J 


z 


ii 


P z 


0 . 0 1 9 1 5 1 5 


RA z 


0.99836604 


J 


z 


3a 


0 z 


0.01 330570 


RA z 


0.99302051 


J 


z 


35 


P z 


0.01859627 


R A z 


0.99760923 


J 


z 


3o 


P z 


0.0! Ra-<693 


R A z 


0.99712106 


J 


z 


37 


P r 


0.0 191 7740 


RA z 


0.99654321 


J 


= 


3d 


P z 


0.01946797 


RA z 


0.9953610 1 


J 


= 


39 


P z 


0.01975='53 


RA z 


0.99505770 


J 


r 


ao 


P r 


0.02004910 


RA z 


0.93(41 i<420 


J 


r 


<11 


P z 


J.02t^ 33967 


R A = 


0.99500880 


J 


= 


42 


p - 


0.')P06S023 


RA z 


0.9917 1690 


J 


= 


as 


P - 


0. 12092 jso 


PA z 


0.9902)067 


J 


- 


au 


P z 


0.021211 S 7 


R A = 


0.94845873 


J 


r 


45 


P z 


0.02150193 


RA z 


0.98642581 


J 


= 


4b 


P z 


0.021 79250 


PA z 


0.98407233 


J 


= 


47 


p z 


0 .02208 50 7 


PA z 


O.P^*! 35125 


J 


T 


4d 


p z 


0.0 22 5 / 3o 3 


PA r 


0.9782222 I 


J 


z 


49 


P z 


p , 0 p p Tb 1 2 ') 


RA z 


0.97162169 
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I 

data output for LINEARfZe SIZE EfpECT UNtAXlAU TENSILE STRESS 
FOR PLATE WITH NOTCH 

WHERE rwe OtSPLACEHEIVT 800NOABY CONDITION IS 0.0100 



half specimen model input 



KC 


= 


0. 


00100000 










5ETA 


s 


0. 


1 ooooooo 










alpha 


z 


10. 


00000000 










PBAR 


z 


0. 


06704500 










N£ 


= 


^2 












ALPHA2 


z 


12. 


OQOOOOOO 










betaa 


z 


0. 


03111025 








J 


z 


1 


P Z 


0.00671700 


RA 


z 


0.99999977 


J 


z 


2 


P Z 


0.00900757 


R<V 


z 


0.99999965 


J 


z 


3 


P z 


0.00P29613 


RA 


z 


0 .99999949 


J 


z 


a 


p = 


0.00956970 


RA 


z 


0 .09599927 


J 


z 


5 


P z 


0.00967927 


RA 


z 


0 . 99999895 


J 


z 


h 


P z 


U.0l0lo9A5 


RA 


z 


0.99999851 


J 


z 


7 


P z 


0.01046040 


RA 


z 


0.99999791 


J 


z 


8 


P z 


0.01075097 


R A 


z 


0.9990971 0 


J 


z 


9 


P z 


0.01 104153 


RA 


z 


0,99999601 


J 


z 


1 0 


P z 


O.Ol 1 53210 


RA 


z 


0.99999454 


J 


z 


1 1 


P z 


O.Ol 1 62267 


RA 


z 


0.99999261 


J 


z 


12 


P z 


0.01 191 323 


RA 


z 


0.99999006 


J 


z 


13 


P z 


0.01220380 


RA 


z 


0.9999067? 


.J 


z 


la 


P z 


0.01249457 


RA 


z 


0.99998239 


J 


z 


15 


P z 


0.01278495 


RA 


z 


0.99997630 


J 


z 


lo 


P z 


0.01 507550 


RA 


z 


0.99996962 


J 


z 


17 


P z 


0.0 1 55660 1 


RA 


z 


0.99996045 


J 


z 


IH 


P z 


0.01 565663 


RA 


z 


0 . 99994860 


J 


z 




P z 


0.01 5^4720 


RA 


z 


0 .9999 5408 


J 


z 


20 


P z 


0.01423777 


PA 


z 


0.99901558 


J 


z 


21 


P z 


0.01432933 


RA 


z 


0.99989242 


J 


z 


32 


P Z 


0.01451890 


RA 


z 


0.99986557 


J 


z 


23 


P z 


0.01510947 


RA 


z 


0.99982777 


J 


z 


2*^ 


P z 


0.0 15*^0005 


PA 


z 


0.99978554 


J 


z 


25 


P z 


0.01569060 


PA 


z 


0.99972912 


J 


z 




p z 


0.015981 1 7 


R A 


z 


0.9996624 1 


J 


z 


27 


P z 


).')lb?7l75 


RA 


z 


0 .99958094 


J 


z 


2® 


P z 


0.01656230 


RA 


z 


0.99948180 


J 


z 


2<^ 


P Z 


0.01685287 


RA 


z 


0.99956158 


J 


z 


30 


p = 


. 0 1 7 1 4 5 1 3 


PA 


z 


0,99921627 


J 


z 


31 


P z 


O.Ol 745400 


RA 


z 


0.99904122 


J 


z 


32 


P z 


0.01772457 


PA 


z 


0.99065100 


J 


z 


33 


P z 


0.01301513 


RA 


z 


0.99857930 


J 


z 


3^ 


P z 


0.01650570 


PA 


z 


0.99027805 


J 


z 


35 


P z 


0.01859627 ■ 


RA 


z 


0.997921 1 7 


J 


z 


36 


P z 


0.01888S83 


RA 


z 


0. 99749601 


J 


z 


37 


P z 


0.0191 7740 


RA 


z 


0,99699403 


J 


z 


36 


P Z 


0.01946797 


R A 


z 


0.99640065 


J 


z 


3^ 


p z 


0.01975853 


RA 


z 


0.99570183 


J 


z 


ao 


P z 


0.02004910 


R A 


z 


0 .99480098 


J 


z 


u 1 


o z 


0.02O S39o7 


RA 


z 


0.99391015 


J 


z 


42 


P z 


0.02065025 


RA 


z 


0.9927940b 


J 


z 


43 


P z 


0.02092060 


RA 


z 


0.991 48300 


J 


z 


44 


P z 


0.021 21 1-57 


PA 


z 


0 . 9aQQi5AS5- 


J 


z 


45 


P z 


0.02150195 


R A 


z 


0.90810624 


I 


z 


46 


P z 


0.021 79250 


RA 


z 


0.9061 3821 


J 


z 


47 


P z 


0.02209 507 


RA 


z 


0.90376971 


J 


z 


48 


P z 


0.02257565 


PA 


z 


0 ,98 1 0 59 4 0 


J 


z 


49 


z 


0.02266420 


Pa 


r 


0 .97789P4S 
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DATA Output fop linearize Size effect uniajiial tensile stress 

FOR PLATE with WOTCH 

WHERE The displacement boundary CONOTTION is 0.01 OO 
ENTIRE specimen MOOEL INPUT 



KC 


= 


0. 


00100000 








beta 


s 


0. 


lOOOOOOO 








Alpha 


= 


to. 


oooooooo 








PSAR 


r 


0. 


08702500 








n£ 


s 


592 










ALPHA2 


: 


12. 


OOOOOOOO 








betaa 


= 


a. 


02887978 








J 


s 


1 


P = 


0.00671 700 


RA z 


0 .99999943 


J 


= 


2 


P = 


0.00900757 


RA z 


0 . 99999015 


J 


: 


3 


P : 


0.00929*^1 5 


RA X 


0.99999676 


J 


: 


a 


P = 


0.00958?70 


RA X 


0.99999621 


J 


s 


5 


P 5 


0.009«7927 


RA X 


0.99999743 


J 


: 


6 


P X 


0.01016965 


RA X 


0.99999636 


J 


: 


7 


P s 


0.01046040 


RA X 


0 . 99909490 


J 






p z 


0.01075097 


PA X 


0.99999292 


J 


= 


Q 


p z 


0.01104153 


RA X 


0.99999025 


J 


: 


10 


p = 


O.Ol 1 3321 0 


RA X 


0.99990666 


J 


s 


1 1 


P s 


O.Ol lb22o7 


RA X 


0.99990195 


J 


s 


12 


P X 


0.01191 323 


RA X 


0.99997572 


J 


= 


13 


P = 


0.01220360 


RA X 


0.99996756 


J 


: 


la 


P X 


0.01249437 


RA X 


0.99995700 


J 


: 


IS 


P X 


0.01276495 


RA X 


0 .99994334 


J 


s 


16 


P X 


0.01 307550 


RA X 


0.99902501 


J 


z 


17 


P X 


0.01 536607 


PA X 


0.9R990342 


J 


z 


Ifl 


p z 


0.0 136S6o3 


RA X 


0.99987496 


J 


z 


1 


P X 


0.013=?4 720 


RA X 


0.99985905 


J 


z 


20 


P X 


0.01423777 


R A X 


0.99979307 


J 


z 


21 


P X 


0.01452633 


RA X 


0.99973733 


J 


z 


22 


P X 


0.01461 690 


PA X 


0.99960689 


J 


z 


23 


P = 


0.01510^47 


PA X 


0.99Q579/JQ 


J 


z 


2'i 


p = 


0.01510OO3 


RA X 


0.99947153 


J 


z 


2S 


p X 


0.0156906'*‘ 


RA X 


0.99935869 


J 


z 


26 


p = 


0.015PH1 1 7 


RA X 


0.9R91 7S07 


J 


z 


?7 


P X 


0.0 1627 1 73 


RA X 


0.99097704 


J 


z 


2H 


P r 


0.01 b5o2‘50 


RA X 


0.99875513 


J 


z 


2^ 


P X 


0.01665267 


R A X 


0.99844! 6 1 


J 


z 


3^» 


P X 


0.01714343 


RA r 


0.99903737 


J 


z 


31 


P X 


0.01743400 


RA X 


0.99766047 


J 


z 


32 


P X 


0.01772457 


RA X 


0.99714793 


J 


z 


33 


P r 


0.01601513 


RA X 


0.99653447 


J 


z 


3^ 


P X 


0.01650570 


R A X 


0.99580245 


J 


z 


35 


P X 


0.01659627 


' RA X 


0.99493150 


J 


z 


36 


P X 


0.016^^66^3 


PA = 


0.99369825 


J 


z 


37 


p X 


0.01917740 


R A X 


0.99267590 


J 


z 


3cH 


P X 


0.01946/97 


R A X 


0.99123383 


J 


z 


39 


p X 


0.01975655 


RA X 


0.9895372 1 


J 


z 


ao 


p X 


0.0200491 0 


R A X 


0.90754645 


J 


z 


ai 


P X 


0.0 20 5 3 9o7 


RA X 


0.P6521678 


J 


z 


a2 


P X 


O.U2065025 


RA X 


0,9«249774 


J 


z 


03 


P X 


0.02092460 


RA X 


0.97933261 


J 


z 


tjy- 


P X 


0.0212M37 


RA X 


0.97565790 


J 


z 


05 


P X 


0.0? 150 193 


PA X 


0.97140315 


J 


z 


U6 


p X 


0.021 79250 


PA X 


0.96648977 


J 


z 


47 


P X 


0 . 02206 30 7 


RA X 


0.96063139 


J 


z 


an 


P X 


0.09237563 


PA X 


9,95433317 


J 


z 


40 


p X 


0.‘»2266 i2o 


‘3 A X 


0.94609176 
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data output for UNEARIZe SIZE EFFECT UNIAXIAL TENSILE STRESS 
FOR Plate without r«oTCH 

WHERE the displacement eoUNOARY CONDITION IS O.OOIO 
entire specimen model input 



KC 




0 . 


001 ooooo 










ft£TA 


= 


0 . 


10000000 










ALPHA 




10 . 


OOOOOOOO 










p6AR 


= 


0 . 


00951200 










N£ 


z 


200 












ALPHA2 


z 


12 . 


OOOOOOOO 










betaa 


z 


0 . 


o«.3>A803 








J 


z 


1 


P = 


0 . 0 n H 7 t 7 ) f) 


RA 


z 


1 .OOOOOOOO 


J 


z 


2 


P 2 


0.0 UP0J757 


PA 


z 


0,99999999 


J 


z 


3 


P 2 


0.:)0929M 1 3 


PA 


z 


0 . 99999999 


J 


z 


a 


^ Z 


0.0093857(1 


9 A 


z 


0,99909999 


J 


z 


5 


P z 


0.00987927 


RA 


z 


0,9999999ft 


J 


z 


6 


P z 


O.OlOlo^ai 


RA 


z 


0 . 9999999 ft 


J 


z 


7 


o z 


0*0104909 0 


R A 


z 


0.99999997 


J 


z 


8 


P z 


0.01075)97 


PA 


z 


0.99999995 


J 


z 


9 


P z 


0.0 1 109 153 


RA 


z 


0 ,99999994 


J 


z 


10 


P z 


0.01133210 


RA 


z 


0.99999991 


J 


z 


1 1 


p z 


O.Ol 1 92297 


R A 


z 


0.99999988 


J 


z 


12 


p = 


0.01191323 


RA 


z 


0.99999084 


J 


z 


13 


p = 


0.01220380 


RA 


z 


0.99999979 


J 


z 


1 


D Z 


0.0 1 2494 37 


RA 


z 


0.99999972 


J 


z 


15 


P z 


0. 0 1278493 


RA 


z 


0.99999963 


J 


z 


1 9 


z 


0.0 1 30755«» 


PA 


z 


0,90909952 


J 


z 


17 


p z 


0.0 1 33990 7 


R A 


z 


0.99999037 


J 


z 


18 


^ z 


O.Ot 395993 


RA 


z 


0 ,669990 1 0 


J 


z 




p r 


0.01394720 


RA 


z 


0.99906896 


J 


z 


20 


P = 


0.01423777 


PA 


z 


0,99999ftb6 


J 


z 


21 


P = 


J. 01452533 


R A 


z 


0.99699830 


J 


z 


22 


P z 


0 . 0 ! a5i 090 


RA 


z 


0.99999784 


J 


z 


25 


o z 


0.01510947 


R4 


z 


0 ,99999 727 


J 


z 


2 '4 


=> = 


0.0 1540 00 3 


0 4 


z 


0.69999657 


J 


z 


25 


p = 


4.41599090 


RA 


z 


0,99909571 


J 


z 


2^ 


P s 


0.01599117 


R A 


z 


0.99999465 


I 


z 


27 


o z 


0.0 1/^271 7 3 


RA 


z 


0.99909336 


J 


z 


2*^ 


P z 


0.0 Io5o230 


PA 


z 


0,999991 70 


J 


z 


29 


t> z 


0. 01985257 


Q A 


z 


0.99998989 


J 


z 


30 


0 z 


0.41 714343 


RA 


z 


0.99998750 


J 


z 


31 


p z 


0.0 1 74 34 0 f) 


R A 


z 


0.99998481 


J 


z 


32 


P z 


0.01 772457 


RA 


z 


0.999981 4ft 


J 


z 


33 


p = 


0.01801513 


RA 


z 


0.99997749 


J 


z 


39 


p = 


0.0 1330574 


R A 


z 


0.90997273 


J 


z 


35 


p = 


0.niP599;?7 


RA 


z 


0.99996705 


J 


z 


39 


0 = 


0.0 1 .58ft 93 3 


RA 


z 


0.99996031 


J 


z 


3 7 


P z 


0.01917740 


R A 


z 


0.99995234 


J 


z 


38 


p - 


n.019'4o797 


R.4 


z 


0.99994291 


J 


z 


39 


p = 


0.01975^53 


RA 


z 


O. 99993 I 8 O 


J 


z 


uo 


p = 


0.02O0491O 


RA 


z 


0.99991875 


J 


z 


ai 


p = 


0.02033297 


RA 


z 


0.99990343 


J 


z 




o r 


0 .02493023 


RA 


z 


0.999«8551 


J 


z 


43 


P = 


O.O20920H4 


RA 


z 


0.906«645R 


J 


z 


aa 


P = 


O.OPl >1 V37 


RA 


z 


0.99954022 


J 


z 


us 


P = 


0 . v) 2 1 5 0 1 9 3 


RA 


z 


0 . 999ft 1 1 89 


J 


z 


49 


P = 


0 .4P 1 79PS0 


R A 


z 


0.99977901 


J 


z 


97 


o = 


4.0220^347 


P A 


z 


0.06674095 


J 


z 


48 


p = 


1 • )2257 3n 3 


P 4 


z 


0 . 99969696 


J 


z 


a<? 


P z 


4 , O?29o420 


C A 


z 


0.09964622 
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Data OOTPUT for linearize size effect UNXAXlAt tensile stress 
FOR Plate without notch 

'WHERE THE OrSPLACEMENT 80UNDART CONDITION IS 0.0100 
ENTIRE SPECIMEN MODEL INPUT 



KC 




0. 


QOtOOOOO 








beta 


S 


0. 


10000000 








alpha 


z 


10. 


OOOOOQOO 








pbar 


z 


0. 


09491000 








ME 


z 


200 










ALPHA2 


z 


12. 


00000000 








betaa 


z 


0. 


04394302 








J 


z 


1 


P = 


0 .0007 1 700 


RA 2 


1 .00000000 


J 


z 


2 


P : 


0.00'^00 7S7 


9A 2 


0 .99999999 


J 


z 


5 


P = 




.?& 2 


0 , 9P9Q9Q9Q 


J 


z 


a 


P 2 


0.')O9S«‘^70 


RA 2 


0.99999999 


J 


z 


5 


P 2 


0.009°79P7 


PA 2 


0.99999990 


J 


z 


b 


P 2 


O.OlO1hP03 


RA = 


0.99999990 


J 


z 


7 


P 2 


0.0 1 0’4bO4O 


RA 2 


0.99999997 


J 


z 


0 


P 2 


0 .(J t 0 7509 7 


PA 2 


0.09999995 


J 


z 


9 


P 2 


o.oiioatss 


RA 2 


0 . 9999999/4 


J 


z 


1 0 


P = 


O.oil 33210 


RA 2 


0.9999P991 


J 


z 


1 1 


P 2 


O.Ol Iba2b7 


RA 2 


0.99999908 


J 


z 


1 2 


P 2 


0.01191523 


RA 2 


0 .99999904 


J 


z 


n 


P 2 


0.01220580 


RA 2 


0.99999979 


J 


z 


1 9 


P 2 


0 . 01 ;? 194 5 7 


RA 2 


0.99999972 


J 


z 


IS 


P 2 


0.01270495 


RA 2 


0.93999905 


J 


z 


1 b 


P 2 


0.01507550 


RA 2 


0.99999952 


J 


z 


t 7 


0 2 


0.01 55bb07 


RA 2 


0.99999957 


J 


z 


\ ^ 


P 2 


0.01 5bSbo3 


RA 2 


0 . 9999991 9 


J 


z 


|9 


P = 


0.01 500720 


RA 2 


0.99999090 


J 


z 


20 


P 2 


}. 01425 7 7 7 


PA 2 


0 .93999000 


1 


z 


21 


P *2 


0.01 452555 


RA = 


0.99999850 


J 


z 


22 


P 2 


0.0 1 18 1^00 


R A 2 


0.99999704 


J 


z 


23 


P 2 


0.01510947 


?A 2 


0.99999727 


.1 


z 


?y 


P 2 


0.01540005 


R A 2 


0 .9Q999057 


J 


z 


2S 


P 2 


0 . ill 5 b 0 0r>4 


R A 2 


0.99909571 


J 


z 


2b 


P - 


O.O 15 O 6 I 17 


R A 2 


0,99999 4o5 


J 


z 


27 


P 2 


O.Ulb271 75 


RA 2 


0.9999953b 


J 


z 


20 


P 2 


0.0to5o25') 


RA 2 


0.999991 79 


J 


z 


2® 


P 2 


O.Olo052‘^7 


RA 2 


0.99990909 


J 


z 


30 


P 2 


^ . 0 1 7 1 4 VI S 


RA 2 


0.9993875^ 


J 


z 


31 


P 2 


O.iM 7 4 54 00 


RA 2 


0.99990401 


J 


z 


32 


P 2 


0.01 7 7 2 4 S 7 


RA 2 


0.99990140 


J 


z 


33 


P 2 


0.01801515 


R A 2 


0.99997749 


J 


z 


30 


P 2 


0.01830570 


R A 2 


0.999R7275 


J 


z 


3S 


P 2 


0.(M85PbP7 


R A 2 


0.9999b7u5 


J 


z 


3b 


P 2 


0 . 0 1 8K8b83 


RA 2 


O.99R90O51 


J 


z 


37 


P 2 


0.0191774) 


RA 2 


0.99995254 


J 


z 


3P 


P 2 


0 . 0 1 9 4 0 7 9 7 


R A s 


0 . 99994291 


J 


z 


JO 


P 2 


0.01975053 


RA 2 


0 . 99993 1 HO 


J 


z 


ao 


P 2 


O.0P00491U 


R A 2 


0.99991075 


J 


z 


u 1 


P 2 


O.OPtj559o7 


RA 2 


0 .99990 54 5 


■J 


z 




P 2 


l.oP(>b 5 02 5 


R A 2 


0.93900551 


J 


z 


u3 


P 2 


0.02092009 


PA 2 


0.99906459 


J 


z 


ua 


P 2 


0. >21 ?1 1 37 


PA 2 


0.99CJH4022 


; 


z 


95 


P 2 


O.OPlSv) 195 


R A 2 


0 . 99981 1 09 


J 


z 


Ub 


P 2 


•).02l 792SO 


RA 2 


0.99977901 


J 


z 


a? 


P 2 


0.O22O8507 


R A 2 


0.9997-4095 


I 


z 


90 


P 2 


0 .022 575n 3 


RA 2 


0 .93909090 


J 


z 


99 


P 2 


O.oPPbo'ipo 


PA 2 


0.93954022 
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APPENDIX C 



COMPUTER PROGRAMS FOR CONVERTING 
FINITE ELEMENT (NIKE2D) PROGRAM OUTPUT TO BE 
STATISTICAL RELIABILITY COMPUTER PROGRAM INPUT 

***************lK'***ifr**^******»***************'****%**«******** 

PROGRAM: CONVERT MTKE20 OUTPUT TO B£ STaTISICAL INPUT 

FOR Ipr. integration entire SPeciMEN /iooel onlv 



************************ ************************************* 



soo 



DIMENSION Xl(59«.) 
BYTE A(30) 

TYP£«, 'ENTER FILE 
READ (5,500) A 
FORMAT 13 OAI ) 
open(unit=s. file 

OPENIUNITsG, file 
WRITE (6, SO) 

N£=S9a 
NP = nE/10-*-1 



,XW(S94), 

CONTAINING rUKE20 

=A, STATUSs'OLO' ) 
s'CC.DAT', STATUS: 



00 100 Msl, 18 

YL(M) = . 1650534 
XW(m):.095 

100 continue 

DO lOl M=19,36 
XL(M)=.0995 
XW(M):.09S 

101 Continue 

00 10? Ms37,l32 
XL(M):. 024075 
XtI(M) = .095 

102 CONTINUE 

00 l03 Mrj33,lS0 
XL(M):.0995 
XW(m):.095 

lOi CONTINUE 

00 104- *1:151 ,160 
XL (Hi:. 1650534 
X«(M)=.09S 

104 continue 

DO 105 M=169,200 
XL(M)=. 1650334 
XW(M):.O2407S 

105 continue 

DO 1 06 M: 20 i , 216 
XL(M 1=.0995 

X\4(M):.02‘f87S 



106 continue 

DO 107 4:217,232 
XL (H)=. 024073 
\N(M):.02487S 

107 continue 

DO 100 m: 233,296 
XL(M):.099S 
XW(M):.02'f875 
100 CONTINUE 

DO 1 09 M=297,316 
XL Cm):. I 658334- 



data OUTPUT* 



MEN* ) 



<W<M) s, 024875 

I O'? continue 

DO no M=317.^28 
XL(M):.»658334 
XW(M> =.024875 

no contimoE 

00 111 «sia9,360 
XLrMi=.099S 
XN(M) =.024875 

1 1 I continue 

DO il2 Ns?61.424 
XI.(M) = .024875 
XX(M) =.024875 

112 continue 

00 113 M=42S,456 
XL (N) =.0995 
XW(m>s.024875 

Hi CONTINUE 

DO 1 14 M=457»4a8 
xl(m)s. 1658334 
X HIM) =.024875 

11 'I continue 

DO llS H=469,576 
XL (M)=. 024875 
XV»(N1 = .0 24875 

115 continue 

Do 118 4=577,592 
XLfM) =.024875 
XN(M' =.024875 

iii> continue 

DO 10 1 = 1 ,NP 
DO 15 7=1,10 

M=10»( 1- ) ) t J ■ 

IF1H.NE.NE) THfN 

IF1J.E0.10) then 

«£AD(5.45) SICY ,SIGZ,SI5X .SrUYZ 

ELSE 

REA0(5,35) SIGX,SIGZ,S1GX,SIGYZ 
END IF 

ELSE 

READ(.S,30) 3IGY,SIGZ, SIGX,SIGYZ 
GO TO 55 

END IF 

WRIT£(6,40) M.XL(M),xW(M),SIGY,SIGZ,SIGX,SIGYZ 

15 CONTINUE 

14 CONTINUE 

SS WRir£(<>,40) M,XL(M) ,XW(M) ,SICY,SIGZ,S1Gy,SIGYZ 

30 FORMAT! l8X,tt£il .0) 

15 FORMAT! I8X,4E1 1 .0.////) 

4') F0RMAT!3X. :S,2F11 .7,4£12.4) 

«S FORMAT! 18< ,4E1 1.0,///'////////) 

50 FORMAT! IX. •ELEMENT NO. ‘ , 4X , • L ' , 1 0X, 'Vl * , 7X, ' S I GY • , 3x , • S I GZ ' 
*,8X, •SIGX* ,7X, 'SIGYZ' ,/) 

STOP 

END 
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►*..»****»*•♦..* ***..4. 

c 

C***«**«*»i» *4 4* 4.. ..*..*. »**.*...****♦* »♦»****** 

c 

C PROGJ^AH: COMVegT NIKE20 OUTPUT 70 &£ STATISTICAL INPUT 
C 

c for 4 PT. integration HALF SPECIMEN MODEL OncY 

C 

Ca****A**A******A****J^A**«*A******A***«***-«***«*««A««**«A*** 

c 

Caaaaaa«aa«aaaaaaa«aaaaaaaaaaa«aaaaaa«aaaaaaaaaaa«aaaa«aa«a« 

c 

OlMCHSION XL('9eS)»’fW(»>28) 

BYTe A(^ 0 ) 

Type., ’ENTER file containimg nike2o ourPUT* 

RFAO(5,500) a 
500 F09M*T( joil ) 

OPEN(UNir=S, filE=a. ST4rus=’0Lt>’ ) 

0PEN(1JNTT = &, FlLc='M.04r’, SFArUSs’NEW’ ) 

wpirE(6,50) 

NE=482 
NEU=ME«4 
HP = MeA./40Yi 

00 t Oo NSI , 100 

XUM) = .0a975 
XW(M)= .028S9 

100 CONTINUE 

DO loi M= 101, 220 
XU f4) = . 024875 
<W(M)=. 028500 

101 continue 

00 102 **s21 1,540 
XL(M'=. 0124575 
XW( fils. 0285 000 

102 CONTTNOe 

00 1^5 Ms541 , 660 

XUN) = . 024875 
XWfuJs. 020500 

105 continue 

00 104 . f1 = 661,76d 
xL(f') = .04975 
xIN(M) = .02050 
lOa CONTINUE 

DO I OS Ms761 , 920 
XV (M)=. 04975 
XN(M)=. 0124575 
105 CO'irCNtJE 



DO 10^ M=92l ,1112 
xL(M>=.O2a075 
XW(M) = .012'j375 
lOb CONfI''it-iE 

00 107 M=1 1 13, 1368 
XUM) = . 0124375 
XWfM)=. 0124375 
107 CONllNUe 

00 l*g H=1369, 1560 
XUM) = . 024875 
XW(MI S.0124375 
KIH continue 

00 109 M: 1561, 1720 
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XL(M) : . 0497S 
XW(M) =.0124375 
109 CONTINUE 

00 no M=1721,1896 
XU(M1=. 0124375 
XW(M)=. 0124375 
no CONTINUE 

00 in M=1897n923 
XL (M)s, 0124375 
XW(M)s. 0124375 

111 continue 

00 10 Isl.NP 
00 15 J»1 . I 0 
00 20 K=l.4- 

N=40 • ( I-l ) +4» ( J-l ) +X 
IF(N.NE.N£4) THEN 
IF(|C.EQ.«) THEN 

iFtj.EQ.to) Then 

R£A0(5,4S) SIGY,3IG7.STGX.SIGY2 
ELSE 

R£AD(5,3S) SIGY.SIGZ.SieX.SIGVT 
EMO IF 
ELSE 

REAO(5.30) SIGT»SIGZ.S1G».SIGTI 
END IF 
ELSE 

REAOfS,30) SICY,S1GZ,SIGX,3TGVZ 
GO TO 55 
ENO IF 

WRITE (6. 40) M»XLCM) rXWCMl.SIGT.STCZ.SIGX.SlGYZ 
20 CONTINUE 

iS continue 

10 CONTINUE 

55 WRITE(6,40) M, XL (M) . XWf N ) , 31 Gr , S IGZ , S I GX , SI GYZ 

30 foRMftT(iax,uen .oi 

35 FORNATI 18X.4E11 .0. /■) 

4Q F0RM4T(3x,rS.2F11 . 7.4E12.4) 

US FORMAT! lex. 4E11 .0, ///////// ) 

50 FORMAT! tX, 'ELEMENT NO , ' . 4X , ' L ' , 1 OX, ' W ' ,7X. ' S I G r • , ax , ' S I GZ ’ 
»,ax. 'SIGX' ,7X, 'SIGYZ',/) 

STOP 

ENO 
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c 

C*******«ft«*-***»*^*****4#^*«*«******#*****«^*«**************** 

c 

c WV3GRAM: CONVERT NIRE21> OUTPUT TO BE STATISTICAL INPUT 

C 

C FOR IPT. INTEGRATION HALF SPECIMEN MODEL ONLt 

C 

C**********************V««****A***^«>*A**A--A*****ft************ 

c 

C*******l^**«********«*ft^»*****#***«»*«***********«****#*1^*A* 

c 

OlrtENSION XL(^ea) FXW(4e2) 

BVT6 A(iO) 

TyPE»,’FIL£ containing M1KE20 OATA OUTPUT' 
ft£AD(5,5O0) A 

500 F0RMaT(30A1) 

OP£N(uMIT=S, fIUE=A. STaTUS='OLO' ) 

. 0PEN(UNIT=6, FJLEs'Sa.OAT* . STATUS* 'NEW ) 

WRITE(&,50) 

NE = if82 
NP=N£/io*1 

00 too M= 1 , 25 
XLfM)=.099S 
Iw(m)=.0570 

100 continue 

DO lOl M=e6,SS 
XL(M)=. 04975 
XW(M)=;oS7 

101 CONTINUE 

00 102 MsSA.IJS 
XUfN) = .O2<v075 
Xw(M'=.OS70 

102 CONTINUE 

00 103 M=136.t65 
XL(M)=.0a97S 
XW(M) = .0570 

103 CONTINUE 

OO 10«. M=166, 190 
XLf M)=.099S 
XW(M)=.057O 

lou continue 

DO lOS M=191 ,230 
XL(M)=.099S 
XV*(M) S.02487S 

105 continue 

DO lOb M=231,279 
XL(N)=. 04975 
yW(M) =.024875 

106 continue 

00 107 M=279,342 
XLtM) =.024875 
XW(M)=. 024875 

107 continue 

00 lOH M=34>J,590 
XL(M) = . 04-975 
XV(fM) = .02A87S 

106 CONTIkuE 

00 109 M=391 ,430 
XL (M) =.0995 
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xwiM)i.0 2a.a75' 

I0<} CONTlHUe 

00 no v»:431,a,74 
XLCMjs.oa'^a^s 
XW(M)=. 024575 
no CONTI Mo£ 

DO 111 M=475,482 
XU(M)=. 024875 
XWfM)=.02487S 
111 COnTINUC 
00 10 Isl.NP 
00 15 J=1 . 10 

Msiom-D+J 

lF(M.Me.N£) THEN 

IF<J.£<J.10) then 

REA0(S,4S) SIGY,3ICI,S1CX,SIGYZ 
ELSE 

REAOtS.SS) 3IG7,SICZ,SIGX,oIGY7 
ENO IF 

ELSE 

REXOCS^oOl SIGY.SIGZ, SIGX,SIGYZ 
60 TO 55 
£N0 IF 

WRITE! 5, 40) M,XL!M) ,XW(M) , SlSf. S IGZ , 3IGX , S IGYZ 
15 CONTINUE 
10 CONTINOe 

55 W(UT£(6.40) M,XL(N),xW(rt),3lGY,SI«2,SIGX,3ICYZ. 

30 F0ftMAT(l8x,4et 1.0) 

35 FORMAT! l8X,4En .0, ////) 
no F0RMAT(3X , I5.2F11 .7r4E12.4) 

45 FORMAT (1 ax, iieil . 0 ,////////////) 

so FORMATfIX, 'ELEMENT NO . ' ,4X , ' L* , I OX, ' W , 7x , • s IG r • , 3x , • S IGZ ' 
•,ax, 'SIGX' ,7X, 'SI6YZ' ,Y) 

STOP 

ENO 
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